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Abstract 

We study from a critical perspective several quantum-electrodynamic phenomena commonly re- 
lated to vacuum electromagnetic (EM) fluctuations in complex media. We compute the resonance 
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shift, the spontaneous emission rate, the local density of states and the van-der-Waals-Casimir 
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pressure in a dielectric medium using a microscopic diagrammatic approach. We find, in agree- 
ment with some recent works, that these effects cannot be attributed to variations on the energy 



of the EM vacuum but to variations of the dielectric self-energy. This energy is the result of the 

00 ' 



interaction of the bare polarizability of the dielectric constituents with the EM fluctuations of an 
actually polarized vacuum. We have found an exact expression for the spectrum of these fluctua- 
tions in a statistically homogeneous dielectric. Those fluctuations turn out to be different to the 
ones of normal radiative modes. It is the latter that carry the zero-point-energy (ZPE). Concern- 
ing spontaneous emission, we clarify the nature of the radiation and the origin of the so-called 
local field factors. Essential discrepancies are found with respect to previous works. We perform 
a detailed analysis of the phenomenon of radiative and non-radiative energy transfer. Analytical 
formulae are given for the decay rate of an interstitial impurity in a Maxwell- Gar nett dielectric and 
for the decay rate of a substantial impurity sited in a large cavity. The construction of the effective 
dielectric constant is found to be a self-consistency problem. The van-der-Waals pressure in a com- 
plex medium is computed in terms of variations of the dielectric self-energy at zero-temperature. 
An additional radiative pressure appears associated to variations of the EM vacuum energy. 
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I. INTRODUCTION 



The concepts of radiative fluctuations and virtual particles are inherent to the nature 
of quantum mechanical processes in the framework of Quantum Field Theory (QFT). The 
actual existence of fluctuations is however questionable as the manner they enter quantum 
phenomena depends on the formalism employed. For instance, one can say that radiative 
corrections renormalize the 'bare' mass of a free electron as long as one can postulate the 
existence of a 'bare' electron in absence of coupling to the electromagnetic (EM) field. Like- 
wise, one can say that zero-point quantum fluctuations of the EM field give rise to a shift 
on the energy of stationary atomic states provided that the actual existence of those states 
could be postulated in the absence of coupling to radiation. However, quantum fields couple 
to each other and such a coupling cannot be switched off. Bare masses and stationary energy 
states are actually unobservable. In addition, if the electromagnetic field appears coupled to 
a stochastic system, eg. a thermal reservoir or a gas of (actual) charged particles randomly 
distributed, additional fluctuations appear and add up to the once previously considered. 
The crucial difference of these with respect to the zero-point fluctuations is that they can 
be modified or even switched off by tuning some external parameter. It is implicit in this 
picture that the spectrum of fluctuations can be expanded as a power series in the free- 
space electromagnetic propagator and the coupling to actual charges, in application of the 
Lippman-Schwinger equation. Many electromagnetic phenomena are explained appealing 
to the actual existence of zero-point fluctuations using quantum or semiclassical approaches 
[46]. That is the case of the Casimir-Polder effect, the Lamb shift, the van der Waals 
forces and the spontaneous decay of excited atoms. Recently [2], it has been pointed out 
that the Casimir effect as formulated in Quantum Electrodynamics (QED) is not actually a 
manifestation of the zero-point vacuum fluctuations of quantum fields but the result of the 
interaction between actual charges and currents -see also [3]. 

In this paper, we are concerned with all the above phenomena in the context of random 
media. We will show that they all can be addressed employing a unified formalism. Our 
study includes the computation of the modified decay rate of an excited emitter and the 
calculation of the dielectric constant of a homogeneous medium. Closely related, we will 
compute the local density of states and the EM energy density of a dielectric. Van-der- 
Waals-Casimir forces and radiative pressure manifest as a response to virtual variations on 
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the total EM energy density due to changes in the spacial distribution of dipoles. To explain 
all these phenomena we will make use of both Classical Optics and QFT formalisms. 

To what Optics concerns, it is a general issue the characterization of a material system 
which couples to radiation by means of its coherent transport properties and the study of the 
decay of unstable local states. With respect to the latter, it is known since the work of Purcell 
[4] that the spontaneous emission rate of an atom, T, in a dielectric host medium depends on 
the interaction of the emitter with the material environment. This is so because the medium 
determines the spectrum of the EM fluctuations which mediate the atom self-interactions 
and hence, its self-energy. In the first place the host medium modifies the density of channels 
into which the atom can radiate -i.e. the Local Density of States (LDOS)- and hence the 
value of T. Second, the integration of the associated self-energy gives rise to a shift in the 
resonance frequency of the atom. Third, the dipole-transition-amplitude gets also modified. 
As a result, the medium is said to renormalize the polarizability of the emitter. 

In addition and complementarily, the medium polarizes the EM vacuum. This reflects in 
the fact that the dispersion relations that EM normal modes satisfy are determined by the 
effective susceptibility of the medium. Complementarity can be seen in that LDOS and V 
depend as well on those parameters which determine the coherent transport features of the 
medium. That is, on the refractive index and the mean free path. As a matter of fact, both 
null transmittance and inhibition of spontaneous emission are expected to occur in photonic 
band gap materials [5]. 

For practical purposes, the understanding of life-times in random media is relevant in the 
context of fluorescence biological imaging [6] and nano-antennas [7] . Also, understanding of 
unconventional coherent transport properties is essential in engineering metamaterials for 
electromagnetic and acoustic waves [8]. 

To what QFT concerns, we will postulate the existence of two distinguishable EM vacua 
attending to the the existence of two different spectra of fluctuations. These are, a source- 
less vacuum in which normal modes propagate and a self-polarization vacuum in which both 
the radiative emission and the photons mediating the self-polarization of the emitter prop- 
agate. While the fluctuations which carry the divergent ZPE live on the former vacuum, 
those which gives rise to the van-der-Waals-Casimir forces live in the latter. In addition, a 
coherent-vacuum in which coherent emission propagates will be identified. 

The paper is organized as follows, 
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In Section II we first analyze the role of electromagnetic fluctuations in the paradigmatic 
quantum process of spontaneous dipole emission. That section serves also to describe the 
features of the framework in which our approach fits. 

Next, we develop our approach in several steps. It is based on a microscopical diagram- 
matic treatment which explodes the seminal works carried out by Foldy [9], Lax [10], Frisch 
[11], Bullough and Hynne [12-14] and Felderhof et al. [15, 16]. Our main goal will be to 
develop a general formalism suitable for treating in the same footing any emission-related 
process in any particular scenario within the framework of linear random media. Thus, the 
next sections are organized in a manner that allows us to address systematically and sepa- 
rately each of the features which characterize a given scenario and a given process. These 
features are, 

A) The nature of the emission. That is, it can be spontaneous emission from an iso- 
lated atom; stimulated emission by an external exciting field on a polarizable particle; a 
combination of both the spontaneous emission of an atom and the emission induced on the 
polarizable molecule which hosts the atom. This matter is addressed in Section III. 

B) The nature of the embedding of the emitter in the host medium. When the emitter is 
itself a host particle or it resides in a host particle we talk of virtual cavity as the medium is 
in fact unaltered by the presence of the emitter. In any other case, the emitter sites within 
a real cavity of radius R. This matter is addressed in Sections IV and V. 

C) The nature of the host medium. For the sake of simplicity, we will restrict ourselves 
to statistically homogeneous and isotropic three dimensional host media. Further classifi- 
cations are made attending to the relation between the wavelength of the radiation, A, the 
correlation length between the host medium constituents, £, and, if it applies, the radius of 
the emitter cavity, R. In the virtual cavity scenario, £ = R and the topology is refereed to 
as cermet topology as the medium is composed by disconnected point dipoles. An effective 
dielectric medium can be defined as long as A > ^. In the latter case the propagation of 
coherent radiation is fully characterized by an effective dielectric constant, e e //. On the 
contrary, any real cavity breaks manifestly translation invariance and no effective medium 
can be defined generally. Nevertheless, if R ^> £, the cavity is said large and the medium is 
seen by the emitter as a continuum. Further, if A ^> £, the medium behaves as an almost- 
effective one w.r.t. the propagation of the coherent radiation emitted from the center of the 
cavity. It is in that case that the topology is that of a simply-connected-non-contractible 
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manifold in which the emitter is placed at the center of a cavity of radius R surrounded by 
a continuous effective medium. All these matters are addressed in Sections V and VI B 3. 

D) The nature of the radiation. In the first place, radiation can be either transverse or 
longitudinal attending to its direction with respect to the propagation wave vector. Second, 
radiation can be coherent and incoherent. The coherent part has two contributions. The 
first one is that emitted directly by the source dipole into the host medium. The second 
one comes from the coherent interference of the latter with the radiation induced in the sur- 
rounding dipoles. Coherent radiation propagates through a coherent- vacuum. Incoherent 
radiation is that which gets dispersed or absorbed. These matters are addressed in Section 
VI together with a study of the process of radiative/non-radiative energy transfer. 

Equipped with the above detailed analysis, we address explicit computations which can be 
performed analytically. In Section VII we concentrate on the virtual cavity scenario. First, 
in Section VII A we compute the spontaneous decay rate of a weakly-polarizable intersti- 
tial impurity within a Maxwell-Garnett dielectric. In Section VII B we compute formally 
the dielectric constant of a Maxwell-Garnett atomic gas for the case that the single-atom- 
polarizabilities are known in free space and non-radiative effects can be ignored. We analyze 
old experiments and propose suitable modifications in order to test our approach. 

We address the real cavity scenario in Section VIII. There, we compute the spontaneous 
decay rate of a weakly-polarizable substantial impurity placed in a large real cavity. 

In Section IX we compare our results with previous ones. Several arguments are given 
to explain why those fail. In particular, we concentrate on explaining why of the erroneous 
use of the so-called local field factors. 

Finally, we explain in Section X the computation of the zero-temperature EM pressure 
in a gas. Both a contribution coming from variations on the matter self-energy and another 
one coming from normal mode fluctuations are found. Only the former is identified with the 
van-der-Waals pressure in agreement with recent interpretations of van-der-Waals-Casimir 
forces. 

We summarize our results in the Conclusions, Section XI. 
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II. THE SPONTANEOUS EMISSION PROCESS IN A DIELECTRIC MEDIUM. 
THE ROLE OF THE ELECTROMAGNETIC VACUUM FLUCTUATIONS 



This section has a double purpose. On the one hand, it offers a comprehensive overview 
of the features which characterize the physical scenarios to which the approach developed in 
the next sections applies. We will clarify to which extend our approach is really microscopical 
and to which extend quantization is taken into account. On the other hand, it describes 
qualitatively some of the subtle points related to the role of the EM vacuum fluctuations. To 
these aims, we analyze qualitatively the phenomenon of spontaneous emission which, since 
the pioneering work by Einstein [17], has been considered paradigmatic in the understanding 
of the EM quantum vacuum. Elements of both Quantum Field Theory and Classical Optics 
will be used. There is little new in this section and our arguments base on well-known results 
of standard textbooks -eg. [18]. Nevertheless, we feel it is necessary in order to appreciate 
the subtle points of the next sections. Reference will be made, with no further explanation 
at this stage, to concepts and results which will appear later. 



A. Fermi's golden rule 

The formula for T is given by Fermi's golden rule. According to it, in second order 
of perturbation theory an excited state of matter localized at f decays via electric dipole 
transition with a rate given by 

I,1IA 

where the dipole approximation has been taken in Eq.(2). \A) is the initial state in which 
the emitter is in an excited atomic stationary state of energy Yiloa with wave function i/ja(x) 
and the electromagnetic field is in its ground state in which there are no actual photons. 
{|A 1lu ia )} is the set of intermediate states in which the emitter is in a lower energy stationary 
state of wave function tpi(x) after the release of a quantum of EM energy, which equals 
huiA = h{ojA — oji) according to the delta function of energy conservation. In free space, 
|7u, JA ) is a transverse photon. (1 is the transition dipole operator and E(r,t) is the electric 
field operator in the Schrodinger picture. Following the introduction by Barnett et al. in 
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[19], by Fourier-transforming the time dependence of E(r,t) into E 0J (r), decomposing the 
latter in terms of creation/annihilation operators and using the completeness of the space 
of EM states, we can eliminate the time dependence in Eq. (2) in favor of uja and arrive at 

r = ^^\M 2 '-^\t I A^-kjm)'- p -, (3) 

where f denotes hermitian conjugate and |f2) s ' p ' is the self-polarization vacuum state which 
is left for a more precise definition. From now on we will simplify matters assuming the 
emitter is a two-level atom -with state labels A, B- with a unique transition dipole matrix 
element JIab = {A\jt\B) = J d 3 x ^a{x)J1{x)'^* b {x). It is important to bear in mind that the 
stationary atomic states so far considered are bound states of actual charges. Those are, 
the atomic nuclei and the atomic electrons. Therefore, they are eigenstates of some effective 
Hamiltonian which already contain electromagnetic interactions as a result of the partial 
integration of 'virtual' photons. No actual photons (i.e. propagating photons) are involved 
in those atomic states and the interaction integrated is essentially electrostatic. Hence, the 
resultant \xab is directly related to the so-called electrostatic polarizability -see Section III B. 
Therefore, once knowing hab the problem becomes one of computing the radiative vacuum 
fluctuations of energy Itwab which couple to the otherwise stationary dipole. 

Applying next the fluctuation-dissipation theorem [20] we can express the electric field 
fluctuations of frequency loja in terms of the imaginary part of a Green's tensor, 

^■({l\k AB (r) % AB {m) s - p - = -h^<Z{Q{r,7-UAB)}. (4) 

Q(f,r;u ab) is the Green's tensor which propagates the electric field from the point dipole 
source which oscillates with frequency uab at r back to itself. That field mediates the in- 
teraction of the dipole with itself. That is, it is the self-polarization field which propagates 
in \Q) s - p \ Alternatively, Q is interpreted as the propagator of virtual photons which are 
emitted and reabsorbed at the dipole source. 

Up to know, we have followed a quite standard procedure. Alternatively, other micro- 
scopical approaches base on generalized optical Bloch equations in the context of second- 
quantization formalism [21-23]. Throughout this paper we will stick to a Green's function 
based formalism. 
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FIG. 1: (a) Feynman rules of the QED field theory. The wavy line denotes the propagator of trans- 
verse (^4_i_) and longitudinal (Ah) photons. The thick dashed line denotes the propagator of charged 
particles. On the right hand side, tree level interaction diagram from Dirac's lagrangian [29]. (b) 
One-loop mass renormalization diagram for charged particles, (c) One-loop vertex renormalization 
diagram, (d) One-loop polarized photon propagator. 

B. Zero-point and in-free-space vacuum fluctuations 

We give here a non-rigorous overview on the derivation of Eqs.(l-3) in free space in terms 
of the interaction between actual charges and EM fluctuations. 

At zero temperature, with no actual charges at all and disregarding virtual ones, photons 
are non-interacting species. They are radiative and hence transverse. Virtual photons are 
depicted as closed wavy lines which stand for frequency modes of the propagator of A± 
with origin and end at the same point. The 'bare' EM vacuum of sourceless modes at zero 

— * — * 

temperature will be denoted by |0). As noted in [24], the propagator of A u equals that for E w 

— * r\ — * 

modulo a prefactor (u/c) . Hence, E = j^A. In application of the fluctuation-dissipation 
theorem we can write, 

(OlJ^J^lO) = ^(0\P ± (Ml\m and (5) 
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\E(f)&(r)\0) = I (0\El(r)E?(r)\0)du 
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J oj 2 Q{G^\r,r,u)}du; = ^J^I J dcouj 3 . (6) 



Only sourceless fluctuations satisfying the transversality condition V • A = contribute in 

— * 

the above equations. Thus, the Coulomb gauge V • A = is a natural choice in this context. 
Longitudinal photons cannot propagate energy as such propagation needs of charged matter 
support. The Green's function in Eq.(6) is that which propagates the electric field of an 
isolated oscillating stationary dipole, 

,2 



V x V x -%l 



G {0 \r,r',uj) = -5^\r-7). (7) 



The left hand side of Eq.(7) is also the Maxwell's equation for the uj-mode of the electric 
field in free space. Although the real part of G^(f, r", uS) diverges for r — > f , its imaginary 
part yields the desired result in Eq.(6) [13]. 

Further, 3{G±\r,r,u)} relates to the LDOS of the electric field which propagate in the 
bare vacuum through [25, 26] 

LD ° ST = -^-J^{Gf\r,f,uj)}}. (8) 

same expression holds for the LDOS of the electromagnetic field A^, LDOS° [24]. It follows 
that the spectrum of fluctuations of E w and LDOS° are proportional, 

LD0S ° = ^{(Ol^^^^l )}' ( 9 ) 
and the total EM energy of the bare vacuum per unit volume reads 

S° = cT 1 J hw LDOS° du 

= || Tr{(0||-(f)lf(f)|0)}d^ = i -^ J duu; 3 . (10) 

The above integral has no cut-off and diverges as ~ u A . That is the so-called zero-point- 
energy (ZPE) and the corresponding fluctuations are referred to as zero-point vacuum fluc- 
tuations. 

When the electromagnetic interaction turns on (i.e. e ^ 0), the theory is that of QED. 
The Feynman rules are the tree- level diagrams in Fig. 1(a). Loop corrections give rise to 
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fermion mass renormalization, interaction vertex renormalization and electric charge renor- 
malization - Figs. 1(6, c, d) respectively. The latter is also responsible for the (virtual) polar- 
ization of the EM vacuum due to charged-field fluctuations. In Fig. 2 we depict a number of 
diagrams contributing to the ZPE of QED. The diagrams in Fig. 2(a) and Fig.2(d) refer to 
the zero-point fluctuations of the non-interacting fields. In the rest, vacuum fluctuations of 
both fields combine. Diagrams (6) and (c) refer to fluctuations in a polarized EM vacuum. 
The fermionic loop of Fig. 1(d) enters the diagrams of the upper row as a one- loop polariza- 
tion function. The pairs of diagrams (6) and (e), (c) and (/) are actually identical. That 
is, from the point of view of the EM field, fermionic loops polarize the EM vacuum. From 
the point of view of charged fields, radiative corrections renormalize their charge. In any 
case, the inclusion of these fluctuations does not make any better the convergence of the 
ZPE. In addition, at finite temperature, the spectrum of thermal fluctuations is Planck's 
and additional photon interaction vertices show up -see [27]. 

Next, let us introduce actual charges. Let us consider an isolated atomic dipole source 




FIG. 2: QED vacuum diagrams contributing to the ZPE. Diagrams (a) and (d) refer to the non- 
interacting theory. The diagrams (b) and (e), (c) and (/) are topologically identical. They give 
rise to a polarization of the EM vacuum as seen by photons (upper row) or to a renormalization 
of the electric charge e 2 as seen by charged fields (lower row) . 

made of a positively charged heavy nucleus and a number of negatively charged electrons in 
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motion around it. In this scenario longitudinal photons are also relevant for they mediate the 
electrostatic interaction between charges. The primitive stationary atomic states incorporate 
actual charges and longitudinal photons. They are eigenstates of the Hamiltonian 

"<> = t^ + Y.-^+t% (") 

i=l i>j | ' i 1 j | i=l I ' 1 1 

where Z is the atomic number; Q is the total charge of the nucleus placed at the center of 
the atom; the subscripts i, j label each electron with renormalized mass and charge m e and 
e respectively; and Pi(xi,t) is the ordinary linear momentum operator of the i th electron. 
The radiative interactions which are not yet integrated in stationary states (i.e. transverse 
photons) enter the (non-relativistic) Hamiltonian of interaction [18], 

~ ~2 

H in t = J2 —A±(r u t) ■ fi(r h t) + -^—-A ± (n,t) ■ A ± (n,t), (12) 

where the magnetic interaction due to spin coupling has been neglected. The dominant 
interaction is given by the first term on the r.h.s. At leading order in e, the representation 
of the self-energy of state A is the second diagram on the l.h.s. of the approx. symbol in 
Fig.5(d 2 ). The dotted line between the points of emission and reabsorbtion there depicts 
the virtual transition state B. For typical frequency transition energies, the wavelength of 
the emitted photons is much longer than the radius of the atoms. In such a case, the first 
term in Eq.(12) reduces to 

—A ± (x i7 t) -pi(xi,t) ~ -ex - E±(xi,t) = -fa ■ E±(xi,t), (13) 

m e c 

where fli is the electrostatic dipole moment operator associated to the i th electron. This is 
the electric dipole approximation, restricted to transverse modes in free space, that is used 
in Eq.(2). This way, the propagator of virtual photons in the self-energy diagrams can be 
reduced to the electric field propagator. That propagator is also the scattering amplitude, 
computed at second-order of perturbation theory in QED, of the process through which an 
excited two-level point dipole transfers a quantum of energy (not necessarily by means of 
radiation) to a non-excited dipole. One can verify for instance in [28] that such a scattering 
amplitude equals the Green's function of Eq.(7), 

r,( Q ) (T} s ^ i-iv-^r^-Bl® (Q\Hi n t\l<?, K )®\A,A)(A,A\® {-y ?tK \ H int \0) <S> \A, B) 



+ 



Hujab — frcq 

(A,B\ (0| H int |7g>) g \B,B)(B,B\ g ( 7g> | H int |0) g \A, B) 

—fkJAB — hcq 

13 



(14) 



where H int is given by Eq.(13), R is the distance between the dipoles and the sum runs over 
all possible values of the momentum q and polarization states k of the only photon in the 
intermediate states, |7g iK ) ® \A, A) and \^fq, K } <8> \B,B). 

In application of the fluctuation-dissipation theorem in-free-space we can write 



t h 



2 ^|^|^{Tr_{Gf(f,r-;^)}}, - * 



f < a, 



3eohe 2 

where stands for the EM self-polarization vacuum seen from the unique dipole in-free- 

space. Because G^\f, r 1 ; ojab) above equals that in Eq.(7) for photons propagating through 
the bare zero-point vacuum, we conclude that the spectrum of fluctuations of forms 
part of the spectrum of fluctuations of |0). 

In a sense, it can be interpreted that zero-point EM fluctuations interact with the dipole 
to make it emit a photon. However, there is no 'borrowing' of vacuum energy as suggested 
in the literature (eg. [28]). The zero-point EM fluctuations and hence the associated ZPE 
remain as depicted in Figs.5(ci, d\). Polarized fluctuations add up on top of those which 
yield the ZPE. On the one hand, fluctuations of the self-polarization vacuum 'interact' with 
the actual dipole making it fluctuate and acquire an additional self-energy. Reciprocally, the 
dipole randomly localized yields an additional effective vertex of interaction which enters 
polarizing the EM fluctuations of the sourceless vacuum In Fourier space, the vertex 

read — (a;/c) 2 V _1 a;o, where a is the electrostatic bare polarizability of the dipole and V the 
total volume of the sample (eventually infinite). Its diagrammatic representation is that in 
Fig. 3 (a). The reason why we depict the bare dipole there as a closed loop is that it is a 
bound state. Momentum and energy can only leak out of it through the coupling to EM 
modes. It behaves w.r.t. the EM fluctuations in as one of the loops of fermionic 

fluctuations. In the process of joining the external fermionic legs of the QED diagrams of 
Fig. 3 (a) in order to get an a loop there is an implicit integration of high frequency degrees 
of freedom. Those are, all the electrostatic interactions which yield the stationary atomic 
state plus high frequency transverse modes. This implies that there is a frequency cut-off 
in all the processes in which a plays the role of an effective vertex. If the radius of the 
dipole is a, the cut-off must be of the order of c/a to preserve consistency with the dipole 
approximation. An analogous reasoning may apply to restrict the interval of frequencies 
which are integrated in the Casimir energy of a system formed by two perfectly conducting 
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parallel plates. In that case, the restriction is imposed to preserve consistency with the 
perfect conductivity approximation [1, 30]. 

So far, we have distinguished between the EM vacuum seen by sourceless photons (i.e. 




A 




± 



o 



FIG. 3: (a) Formal integration of electrostatic and rest of high frequency modes in the stationary 
atomic state A and bare electrostatic polarizability cto- (b) Effective interaction vertex in the 
dipole approximation after integration of high frequency modes, (c) Vacuum polarization diagram 
in presence of a unique dipole. 

normal EM modes), |^)f^, and that seen by the photons whose source is the dipole, 
This distinction is based on the double role played by loops of actual charged-fields and 
loops of virtual photons. Those are, the renormalization of the bare polarizability and the 
polarization of the EM vacuum. The same as shown in Fig. 2 for the ZP QED vacuum, both 
effects are complementary. The EM vacuum seen from the emitter was identified above, 
at least partially, with the bare zero-point vacuum. That is, the radiation emitted by the 
dipole sees empty space in its way. Likewise, self-polarization photons see empty space from 
and back to the dipole. However, the restrictions uo < c/a must apply in order to consider 
the emitter as a point dipole. Hence, the equivalence r — r 1 implicit in Eqs.(5-10) has been 
substituted by the limit \f — f*\ C a in Eq.(15) to set explicitly the limitation in spatial 
resolution. Therefore we will write = |0)| u ,< c / a . The associated LDOS is 



In an analogous manner that fermionic loops renormalize the fine- structure constant ~ e 2 
in QED through its coupling to radiation, the ao-loop renormalizes the value of «o itself. 
This self-polarization effect is diagrammatically given by the geometric series in Fig. 6. It 



LDOS' 



emission 



ID 



4ll 2 C 2 w<c/a 



(15) 
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yields the renormalized polarizability a in free space. Note that the photon loops in Fig. 6(c) 
denote Q u3 {f ) r) which equals G^°\f,f;uj) in free space, where r is the position vector of the 
dipole. 

On the contrary, the vacuum that sourceless normal modes see, |0)f ^5, is polarized by 
the presence of the dipole. In other words, radiation traveling from infinity and normal 
modes can meet a randomly distributed scatterer in their way. Its EM fluctuations are 
made of two additive contributions. The first one yields the divergent ZPE. The second one 
is depicted by those diagrams in which virtual photons of frequency less than c/a scatter 
with the dipole of 'dressed' polarizability a -see Fig. 4. Note the difference with respect to 
those in Fig. 6(c). In the diagrams of the polarized sourceless EM vacuum two propagators 
attach to the dipole from below. That is so because the whole series in Fig. 4 denote 
the bulk loop-propagator, Q{G±(r,r)}, where r is any point. In contrast, the loops in 
Fig. 6(c) stand for ^{^(r, r)} = ^s{G^(f,r)} for f being the position vector of the dipole. 
The physical distinction between G and Q will be explained in due course. The LDOS for 
sourceless-normal modes reads 

(16) 



LDOS 1D - Lp QS^+ 2V 47r2c2 



ui<c/a 

Note that, because a depends on the fluctuations of LDOS^™ sslora is implicitly ac- 
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FIG. 4: Diagrammatical representation of LDOSl°^ rceless . The two terms on the r.h.s. correspond 
to the integrands of the two terms on the r.h.s. of Eq.(16) respectively. The geometrical series of 
the second term stand for the states multiply-polarized by the presence of the unique dipole. 

counted for in the second term on the r.h.s of Eq.(16). 

At leading order, the counterpart of the addition of the diagrams of Figs. 2 (a), (b) of the 
virtually polarized vacuum is the sum on the r.h.s of the approx. symbol in Fig.5(cq) of 
the actually polarized vacuum. Likewise, the counterpart of the addition of the diagrams 
in Figs.2(<i), (e) for the renormalization of the fine structure constant, is the sum of the 
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(c t ) Microscopic (QHD) representation of the combination process as seen from sourceless photons 
(c.%) Microscopic (QED) representation of the combination process as seen from in-free-space the dipole 




(</ ) Leading order vacuum polarization diagrams for sourceless slates - EM vacuum energy 
(d 2 ) Polarizability renorma ligation at first order in radiative corrections ~ dipole self-energy 

FIG. 5: Interpretation of the 'recombination' process between a point dipole and zero-point EM 
vacuum fluctuations, (a) Zero-point vacuum fluctuation loop of the electromagnetic field. The 
equivalence in Eq.(5) has been depicted. (b) Feynman diagrams contributing to a stationary atomic 
state and its equivalence after integration of high frequency modes as in Fig. 3. (c±) Microscopic 
representation of the recombination process leading to polarized EM vacuum fluctuations in the 
presence of an actual dipole. (02) Microscopic representation of the recombination process leading 
to an actual fluctuating dipole. Equivalent interpretation: renormalization of the dipole polariz- 
ability through radiative corrections. In (^1,2), same phenomena as in (dip) using the effective 
nomenclature of Fig3. 

diagrams on the r.h.s of the approx. symbol in Fig.5(<i 2 ) for the renormalization of the 
polarizability. However, the equivalence between the last diagram on Fig.5(cq) and the last 
diagram on Fig.5(d2) does not hold because the 'fermionic loop' in Figs. 2(d), (e) turn into 
an actual dipole in Fig. 5. Although both sourceless and sourced photons see statistically 
homogeneous vacua, \Q) s - p - is attached to actual dipoles while \Q) Smlm is translation invariant. 
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As a matter of fact, the series in Fig.5(d 2 ) stands for the self-energy of the actual dipole. 

We next overview the renormalization of the polarizability in-free-space. As mentioned 
above, the imaginary part of (r, r) yields the right result consistent with the quantum 
radiative corrections which dress up the electrostatic polarizability [13]. However, it contains 
both longitudinal and real transverse parts which diverge. The phenomenological regular- 
ization scheme that we use in Section III B follows that of [35]. It allows us to attribute 
a physical interpretation to the divergences. The longitudinal divergence is associated to 
the longitudinal modes which are integrated out in the original stationary state. Therefore, 
it relates to /iab and ojab- The associated electrostatic polarizability ao is proportional to 
fJ'AB ^n a two-level atom. The real transverse divergence is interpreted in terms of a shift 
in the transition energy with respect to the stationary original state. That gives rise to a 

2 

static shifted polarizability, oP stat (u) = a i , «o being real. The shift in energy is the 
analog to the Lamb shift in the Hydrogen atom. That scheme of renormalization yields the 

2 

Lorentzian-type polarizability in free space a{uo) ~ ^ 2 a °^2t[ iro • There have been proposed 
however several schemes to accomplish the computation of a in-free-space for idealized two 
and three- level atoms. We refer to the reader to the works of Berman, Boyd and Milonni [31] 
for a computation based on Heisenberg and Schrddinger's pictures; to that of Barnett and 
Loudon [32] for a Green function approach; and to that of Bialynicki-Birula and Sowinski for 
an approach based on Feynman diagrams [33]. Even in the latter work, the authors begin 
with a priori modeled electronic bound states and perform a dimensional reduction from 
three to zero spatial dimensions. It has been pointed out in [31] that a Lorentzian profile 
is an approximation not fully justified and it has been found in [34] that other more com- 
plicated form can be more suitable. We will not attempt here to analyze those approaches 
and we will stick to the phenomenological Lorentzian profile. Also, we will adapt in Section 
III a phenomenological treatment appealing to the existence of a classical effective dielec- 
tric constant within the scatterers. That procedure is in all equivalent to the actual QED 
treatment once divergences in (r , r , to) are conveniently regularized. 

C. In-random- medium vacuum fluctuations 

In a random medium, on top of the in-free-space fluctuations, additional fluctuations are 
induced by the multiple scattering of light with the host scatterers and the stochasticity of 
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FIG. 6: Diagrammatical representation of the geometrical series which give rise to the radiative 
renormalization of the electrostatic polarizability of a two-level atom in free space. In (a, b) it is 
written as a succession of coupled virtual states. In (c), the effective rules in Fig.3 are used. In 
(b, c), the dipole approximation is implicit. 

the medium itself. The additional fluctuations affect both the polarization of the sourceless 
EM vacuum and the self-energy of matter. 

The ondulatory nature of light implies that fluctuations appear as considering the in- 
terference between all the different paths that a virtual photon can follow in its way. The 
particularity of a random medium with respect to the free space is that the different paths 
are not statistically equivalent as a result of the spatial correlations between host scatterers. 
Therefore, each path has a statistical weight and fluctuations are affected by the randomness 
of the ensemble of scatterer configurations. In turn, this gives rise to a stochastic equation 
whose Green function incorporates the EM fluctuations we seek for. Purely thermal fluc- 
tuations of the electromagnetic field will not be considered in this paper. However, it can 
be intended that the fluctuations induced by spatial correlations have an indirect thermal 
origin. That is, it is the classical degrees of freedom which characterize the configuration 
of host scatterers, vector positions, velocities and angular momenta, that obey a thermal 
distribution. In our approach, atomic orbitals do not overlap and the separation between the 
emitter and the host scatterers and the separation between host scatterers themselves are 
determined by the range and strength of the cohesive forces -but for the dipole interaction- 
as well as the temperature. In fluids such atomic gases, one can think of that force as the 
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short-ranged repulsive force between the external electrons of the atoms. In experiments, 
the range of this force can be tuned by ionazing the atoms or molecules [71]. For practical 
uses, throughout this paper we will consider an effective rigid exclusion volume around both 
the emitter and the host scatterers. In the construction of the dielectric constant in Section 
VII B we will only assume knowledge of the in-free-space polarizabilities of the atomic spaces 
considered as isolated. In crystals and tight-binding models in general, atomic orbitals nat- 
urally overlap and the electronic band structure of the system is completely different to 
that of the atomic species. The assignment of a 'new' single polarizability to each atomic 
site might be still possible provided that the only remaining interaction between individual 
atoms is effectively dipole-like. Nevertheless, for the case that the emitter is a distinguish- 
able dipole, the only relevant information for studying the influence of the medium in its 
decay rate is that of the spatial embedding of the emitter into the medium provided the 
medium is correctly described by an effective dielectric constant -see Section IV B. 

We summarize here the basis of our approach. In the process of dressing up the elec- 
trostatic polarizability ao, all the virtual photons running inside each dipole are assumed 
integrated out. The interaction between atoms is dipole-like provided atomic orbitals do not 
overlap. The 'bare' propagator of such an interaction is the one given by Eq.(14) and the 
'bare' interaction vertex is given by the strength coupling — k^pao. This way, in-free-space 
radiative fluctuations and those induced by the host medium can be incorporated in the same 
footing. Quantization enters at the level of atomic polarizabilities, dipole-dipole interaction 
and EM fluctuations. In our approach, matter fluctuations enter classically through the 
implementation of the correlation functions between scatterers. These correlation functions 
are affected neither by the dipole-dipole nor by the radiation-dipole interaction. Thus, the 
coupling between matter and radiation is assumed weak unless specified otherwise. Explicit 
treatment of quantum long wavelength excitons, typical of tight-binding models, is disre- 
garded. The reader is referred to the pioneer works by Fano and Hopfield [36, 37] for more 
convenient approaches to highly ordered systems. As long as collective quantum excitons 
can be neglected, and in the weak-coupling regime between medium and radiation, there 
is no need to develop any further macroscopic quantization scheme. In particular, there is 
no need to postulate any additional noise polarization operator -compare with [38]. Only 
for the case that the emitter is a substantial impurity we will assume any other fluctuation 
incorporated in the given dielectric function of the host medium without further explicit 
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computations -see Sections V,VIII. Additional non-radiative effects, such as collisional ef- 
fects, can be incorporated through convenient renormalization of the vertex — k^pao. 

Next, let us consider formally the host medium as an stochastic configuration of dipoles, 
where stochasticity concerns the classical variables which determine the configuration. In 
order to treat the EM fluctuations induced by the classical matter fluctuations in the same 
footing as the rest, it is necessary for them to fulfill some conditions. The Green's func- 
tion which appears in the fluctuation-dissipation theorem of Eq.(4) refers to an effective 
and stationary configuration. That implies conditions over the dynamics of the ensemble of 
scatterers with respect to that of the emission process and photon propagation. In the first 
place the host medium must pass through all the accessible configurations in a time scale 
much less than T -1 . That is, the relaxation time of the ensemble, r, satisfies r <C T^ 1 . This 
guarantees that the degrees of freedom of the host medium are the fast variables which can 
be integrated out. Second, stationarity demands that r 3> u^ 1 . In addition, it is implicit in 
the derivation of Fermi's Golden rule that 3> lu^ 1 , which is a weak-coupling condition 
on the coupling radiation-emitter. That way, the spontaneous emission process in such a 
random medium adjusts to the Born-Markov approximation [39]. We can also appeal to 
the ergodic character of the electric field so that ensemble averages are equivalent to time 
averages [40, 41]. 

We follow next a similar treatment to that of [41] in the introduction of the matrix den- 
sity operator. \£l) s - p - is not |0)| aJ<c / a any more as the surrounding dipoles polarize |f2) s - p - 
too. There exists a one-to-one correspondence between the spatial configurations of scatter- 
ers and the corresponding polarized vacuum states. Let us denote those vacuum states by 
{\(pm) s ' p '}- The density matrix operator M is a diagonal matrix whose entries M mn are the 
statistical weights of the vacuum states, 

M EE Ys M ™n\K) S - P - S - P -{<t>n\, M mn OC 5 mn . (17) 

m,n 

In the simplest case, M mn is a Boltzmann weight. Application of the fluctuation-dissipation 
theorem over each EM vacuum yields 

s - p i<P m \ME^{r)E^\r)\<p n ) s - p - = f; cu AB )}M mn , (18) 

where the Green function <5 m „(f, r; ujab) corresponds to that of the self-polarization field of 
a point dipole sited at r surrounded by the m th configuration of scatterers. As the medium 
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is stochastic and statistically homogeneous and isotropic, the final expression for the decay 
rate so obtained reads 

T = -^^\»AB\ 2 %{Tr{M mn ® mn (f,r;u; AB )}}, (19) 

where Tr denotes the trace operator. The action of Tr must be intended not only over the 
spatial indices of the tensor propagator, but also over those indices corresponding to the 
space of EM vacuum configurations. The EM vacuum \£l) s - p - is this way a superposition of 
states. Following [41], each of those vacuum states \4> m ) s ' p ' is weighted by the root-squared 
of the coefficients of the density matrix, \il) s - p - = J2 m v / ^mm|0m) s ' p '- In terms of it, we can 
write 

s - p i<p m \ME^ AB (r)E^ AB \r)\<p n ) s - p - = s - p - (Q\E UAB ' (r)E WAB \r)\Q) s - p - , (20) 

m,n 

Summing over the indices denoting the vacuum space configurations we obtain 

G(r, f; uab) = ^2 & mn (r, f; u A B)M mn . (21) 

m,n 

Note that in the above development we are implicitly employing a picture in which the 
electric field operator remains unperturbed under changes in the host medium while it is 
the vacuum state that changes. 

Alternatively, Q{r, t^uab) is the Green function of the stochastic equation for the self- 
polarization field of a point dipole surrounded by a stochastic configuration of host scatterers. 
In turn, Q is determined by a stochastic kernel that relates to the electric susceptibility 
tensor, x, which characterizes the medium. Finally, we can write jl • Q(r,r;uiAB) • ft = 
||/i| 2 Tr_{^(r, r; uab)}, where now the trace operator Tr_ only acts over spatial indices and 
the factor 1/3 amounts to averaging over the three equivalent directions so that 



2uj ab 



- 3£o ^-I 2 ^-{S>>>.b)}}, («) 
and correspondingly, LDOS™ io " = --^3{Tr_{^(f,f;cj)}). (23) 

Z71C I J 

In summary, during a time of the order of T -1 , self-polarizing photons of frequency ujab 
explore all the accessible states compatible with both energy conservation requirements and 
environmental constraints. To each configuration of the environment corresponds a different 
set of possible photonic paths. They all add up, interfere, and give rise to a spectrum of 
fluctuations ~ $5|Tr_{^(r, r; c<j)}| which obey certain stochastic equation. 
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The propagator Q(r,r;ui) must not be confused with the Green's function of the normal 
modes which propagate in The latter is denoted throughout this paper by G and 

is referred to as Dyson or bulk propagator. The LDOS associated to G is LDOS^ ource/ess . 
These matters will be addressed in detail in the next Section. 



D. Longitudinal vs. transverse emission 



As mentioned in Subsection II B, it is only in free-space that Q(f, r; a;) equals the propaga- 
tor of Maxwell's equation for the electric field in free space, G^(f, f; uo). G^(f; uo) consists 
of an electrostatic (Coulombian) dipole field propagator, 



stat 



(r) 



4n r 



(24) 



plus a radiation field propagator, 
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(25) 



-Anr /eg —Anr 
with k = cu/c. In reciprocal space and for isotropic systems, any tensor can be decomposed 
into a longitudinal and a transverse part with respect to the propagation direction of the 
wave vector k. That is, Q(k) = G±(k)(t — k®k) + Q\\ (k)k <g> k, where k is the unitary vector 
parallel to k. In free space, 

\GZXr).(t-k®k) + G?lXr).k®k] 



&°\k)= J d 3 re ir G(°»(r) = Jd 



3 r e ik-r 



Gf\k)(I-k®k)+G\l' > (k)k®k, 



(26) 



with 



Gf(*) 
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h 2 — h 2 ' 



(27) 



While the radiative component is fully transverse, the electrostatic one is fully longitudinal. 
In most of this work we will deal with the Fourier transform of G^°\ The reason being that 
it is easier to keep track of the longitudinal and transverse contributions to the emission by 
direct computation of 



0Ur,r) = iTr_{^(f,f)}I=V 
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d s k 



<3\\{k) 



(28) 
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where the factor 2 in front of G±(k) comes from Tr_{I — k <8> k} and stands for the two 
transverse modes. 

In a generic medium, there does not exist direct identification between radiative and 
transverse emission and neither does it between electrostatic and longitudinal emission, but 
for in free space. In free space the whole radiation consists of transverse modes which are all 
coherent propagating modes as they satisfy the free-space (on-shell) dispersion relation k 2 = 
k 2 .. That is, they are poles of G^(fc). Because G^ does not contain poles and longitudinal 
modes do not couple to transverse radiative ones in free space, the above identifications 
are possible. On the contrary, the presence of a polarizable environment gives rise to non- 
coherent radiation as well as absorbtion which contain both longitudinal and transverse 
modes. Unless a very specific configuration of scatterers holds -eg. a dipole chain [42]- the 
coherent propagation is entirely transverse and corresponds to the far-field coherent signal 
which would be received by a distant antenna in the medium. The incoherent radiation is the 
noise signal received by that antenna. It has its origin in the coupling of 'bare' radiative- 
transverse modes to 'bare' electrostatic-longitudinal modes and contains both transverse 
and longitudinal effective modes. The coupling takes place both at the host dipoles and 
the cavity surface. Note also that, differently to the self-polarization field in free space that 
appears in Eq.(15) which is purely transverse, the one in Eq.(20) and thereafter contains 
both transverse and longitudinal modes. The reason being that while all the longitudinal 
modes have been integrated out in the bare polarizability of the unique dipole in free space, 
still long wavelength longitudinal modes mediate the interaction (i.e. the induction) between 
distant dipoles in a random medium. Nevertheless, the fact that the emission be longitudinal 
does not imply that it is not radiative as erroneously interpreted in some works. 

In free space, 



3 



{Tr_{^(f,f)}}| /ree = 3{Tr_{G^>1}} 

= ™{J^ksh?} = -*><*"> (29) 

the emission is entirely radiative and, in application of Eq.(15), we can write T = 

3 

3^0 fic 3 I^abI 2 , where the script refers to the use of in-free-space parameters. The real parts 
of both 7^°' ) = J J^i2 G ( f > (k) and 7jj°' ) = J j^sG\\(k) diverge and are needed of regulariza- 
tion. Because they are the zero-order terms of the perturbative expansion of Tr{Q w (r, r)}, 
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we can write 

^mr,r» = J ^2G±(k) + J ^G\\(k) (30) 
ee 2 7 ™- + 7 ™- = [2 7 f + 7 J 0) ] + 2 7± + 7|| , (31) 

where 2 7j _ and 7 y are the divergenceless pieces. Physically, the 7 -factors account for the 
dipole self-energy associated to intermediate EM states of frequency uo. The self-polarization 
field propagated with Q w is also referred to as local field. The imaginary part of the 7 -factors 
relates to LDOS* mission through 

LDO g~n = _^_3 {2(7 (0) + 7±) + 7||} _ (32) 

It is also possible to relate T with the power emitted in the decay process. That is, 
because such power is due to the self-interaction of the dipole with its own field, we can 
write 

W, = ^p{/2- E* exc } = ^%{fi-g: r jr-f) ■ ff} 



CO 3 
res 



/i| 2 ^{Tr{^ res (f,f)}}, (33) 



6c 2 eo 

where in the second equality the exciting field E* xc has been identified with the self- 

2 _ 

polarization field, ^^G^ res if, r) • ft*- The values of uj res and jl are those of the resonance 
frequency and the transition dipole in the medium (which may differ from the in-free-space 
values, ujab, A* as)- Comparing Eq.(23) with Eq.(33) we can write T M = jj^WV 



III. SPONTANEOUS VS. STIMULATED EMISSION 



A. Stimulated emission, W w 

Let us considered first that the emitter is a dipole stimulated by an stationary external 
field. The exciting field, E exc (r) = Eq(t), oscillates in time with frequency uo far from any 
resonance frequency of the dipole. In this case, the processes of absorbtion and emission 
of radiation by the emitter become stationary after a relaxation time of the order of T -1 , 
which is supposed much greater than the relaxation time of the host medium. Therefore, 
the process through which the absorbtion and emission become stationary is Markovian 
with respect to the dynamics of the host medium. The stochastic computation of Q is so 
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justified as in the spontaneous emission process. In order to regularize the divergences in 
we consider the dipole as a spherical scatterer of radius a and relative dielectric constant 
e e . This model implies a classical regularization scheme. Again, the dipole approximation 
requires a <C k$ 1 , with k = u/c, A = 2n/k Q . Formally, the emitted power reads, 

W u = J d 3 r 0(r - a)p-(f) • E%*(r)}, (34) 

where fF{r) is the density of dipole moment induced, which is proportional to Eft (r) in 
our linear and small particle approximation and is affected by self-polarization effects. We 
introduce the self-polarizing field through the insertion of appropriate Green's functions in 
the above expression, 

p*{f) = J d 3 r"0(r - a)e Xej dVG^f, 7) ■ [G ( %\f, f") ■ E%(7'), (35) 

W u = |s{y d 3 r Xe Q(r-a) J dVd 3 r"G w (f , f) 

• [G^]~ 1 (f',r") ■ Eq{t") ■ i?o*(r) j. (36) 

In these expressions, Xe = (e e — 1) is the relative electrostatic susceptibility of the emitter 
-not to be confused with the susceptibility of the random medium- and 

oo „ 

G w (rO » G(°)(f) ^ [-A; 2 Xe / 0( V - a)£»d 3 V ] m (37) 

m=0 

is the propagator which makes account of the infinite number of self-polarization cycles which 
give rise to radiative corrections. G U! (f—r > ) propagates virtual photons from a point r 1 inside 
the emitter back to another point r also within the emitter. All the equations above become 
simple in the small particle limit, a <C kg 1 , for the electric field is nearly uniform within the 
emitter and so are the density of dipole moment and the propagator G u (f, f). The n-points 
irreducible diagrams which enter the computation of G u (f) can be approximated by the 
series of Fig. 7(6) in which the two-point correlation functions 9(r — a) appear consecutively 
as factors of a product. That way the corresponding integrals appear disentangled and the 
corresponding series becomes geometrical. The underlying approximation is J d 3 r 0(r — 



a)g u (r) ~ ^a 3 ^(0) = ^a 3 2^' + 7j Uj + 2-f ± + 7y |I We already mentioned that 27 
and 7||°' ) are divergent. Those divergences are cured in our classical model by the presence 
of the finite radius a. Since the limit of the integral Lim{ f d 3 r 0(r — a)G^ at (r)} = -^1 as 
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koa — > is conditionally convergent, it is the Heviside function 0(r — a) of the integrand 
which we used to model the spherical shape of the dipole which yields the finite value -^1 
[43]. Any other geometry would give a different numerical value. By equating that result 
with ^a 3 7,{ 0) p (leaving $t{2^ ] } still free) we obtain 7 j 0) = (^a 3 ^) -1 . The net effect 
of this regularization procedure is the dressing up of the single particle susceptibility with 
all the in-free-space electrostatic corrections. This procedure is depicted in Fig.7(c). Its 
quantum counterpart is the integration of the electrostatic interactions which give rise to 
an atomic bound state in a two- level atom. That way we can define Xe = 7^2 Xe and obtain 
the bare electrostatic polarizability a = 47ra 3 ^i. We emphasize that it is the electrostatic 
polarizability «o which really has physical meaning regardless of the regularization scheme 
applied in its computation. We refer to section III of [35] for a comprehensive summary 
of regularization methods. With the above definitions we can rewrite Eq.(36) in terms of 



♦ s — + C^i + £X^+^++^ + -= — [ 1 + i/3Tr£^i}]" 
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FIG. 7: (a) Feynman's rules for the classical regularization scheme of Section III A. (b) Dia- 
grammatic representation of Eq.(37). (c) Diagrammatic representation of the dressing up of Xe 
leading to ceo- Approximation symbols denote that the field within the emitter is taken uniform. 
(d) Diagrammatic representation of a the action of a self-polarization cycle on the electrostatic 
polarizability. 

electrostatically renormalized operators, 

W u = ^T%{f d 3 r-XeQ(r-a) J dVdV'G^r, r) 

■ [GWJ-V ,f") • E%{?') ■ £ "*(r)}, (38) 
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where 

oo 

G w (r,0 = G< >(r,O ^(-A; 2 a ) m 3- m (2 7 | 0) +2 7± + 7! |) m . 

m=0 

Note that longitudinal and transverse modes couple to each other in the series of Eq.(39) 
as scattering takes place off the dipole surface. This is the classic phenomenological analog 
to the quantum coupling of longitudinal and transverse modes of the diagrams of Fig. 5 (02). 
Finally, in function of the 7-factors the power emitted and absorbed by the induced dipole 
reads 



w = ^ r 



2 ~ll + I]fcg Qo [2 7 f +27x + 7||] J 

{ ^ Q {2l lot. + Tot. } (4Q) 

^|l + iHanf27° ) + 27 l + 7„ll 2 " 



6c 2 



3 S{a } 



■}\E%\ 2 . (41) 



^ 2 |l + |A; 2 ao[27f +2 7 ± + 7||]P 
The term in Eq.(41) corresponds to the power absorbed within the emitter. The term in 
Eq.(40) corresponds to the power radiated into the medium. The former is non-zero only 
if e e contains an imaginary part [96]. We can write Eq.(39) in terms of a renormalized 
polarizability a as W w = ^3{<x% • (Eft)*} = ^|££| 2 3{<5} with 

« = mT^ > ( 42 ) 

1 + |/c 2 a 7r + |^o«o[2 7 ± + 7||] 



so that Eqs. (40,41) can be written also as 



^" 6c 2 1 1 



3 %{a } 



c J{27 - +7r . } __ ! ^|. (43) 



B. Relation between induced and spontaneous emission for the Lorentz dipole 
model, F a 

The expression for a still contains a divergent term in the denominator, |/cqO:o7j ' ) - While 
its imaginary part is convergent and gives rise to the usual radiative corrections for an off- 
resonance point dipole in vacuum, a ~ a [l — -^k^a^ 1 , its real part diverges and needs of 
regularization. Following the regularization scheme of [35] we can take advantage of this and 
incorporate the divergence into a resonance in a. The procedure must be compatible with the 
regularization scheme used for a and it is our choice that it reproduces the phenomenological 
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Lorentzian model (L) for the polarizability of a two-level atom, 

a L (k) = ^ok 2 r JkL -~k 2 - ir a ~k 3 /(ck 2 res )]-\ (44) 

where k res is the resonance wave number and 5o is the 'renormalized' value of the electrostatic 
polarizability of the emitter, both evaluated within the medium. We proceed by equating 
Eq.(42) and Eq.(44) in absence of absorbtion within the emitter and with all the parameters 
evaluated in free-space - which is denoted with script 0, 

a [l - ^a k 3 + ^ PK{2 7 | 0) }]- 1 = a„[l - k 2 /k 2 - iT^]- 1 . (45) 

From Eq.(45) we identify U{2^} = and T = ca ^o/67r. We refer to a° stat = a [l — 
k 2 jk\\~ x as the 'shifted' electrostatic polarizability at frequency u = ck. By consistency 
with the decay rate of a dipole with transition amplitude /xo in free-space according to 
Fermi's Golden rule, Tn = „ k ° k \uo\ 2 , we find the relation an = 2 ^°j . The latter can 
be found in textbooks (eg. [44, 45]). It is important to emphasize that, although the 
parametrization of a in Eq.(44) is phenomenological and questionable on the basis of some 
two- level atom quantum models [34], the renormalization procedure after regularization of 
divergences is independent of such parametrization. Our renormalization scheme bases on 
the renormalization of the photon propagator in which, after regularization, —k 2 pa° stat enters 
as a point vertex. Hence, although the computed values of all the renormalized parameters 
depend on the functional form attributed to the in-free-space polarizability, the local density 
of states computed through Eq.(23) is independent of such arbitrary choice provided that 
the value of Oi° stat is the correct one. 

In the host medium, the values of T, k res and oq get renormalized according to the 
equations, 

r« = -^^^{27r+7r}i^ res 

= -r |^{2 7 r + 7 p-}l^ es , (46) 

where k res is a real non-negative root of the equation 

(k/ko) 2 ~ 1 = \^k 2 ^{2 1± + 7 ||}k =fcres , (47) 
and oq = a (k /k res ) 2 . (48) 
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Eq.(47) is the Lamb-shift of the resonance frequency, which is as a result of the variation 
of the real part of the dipole self-energy [47]. In addition, Eqs. (46-48) together with the 
requirement of consistency with Fermi's Golden rule imply that /i gets renormalized with 
respect to /i . 



C. Combination of spontaneous and induced emission, r° 

Finally, consider the spontaneous emission of a point dipole like that in Eq.(22), but 
now with an additional bare electrostatic polarizability a . The situation is analogous to 
that of a fluorescent atom with transition dipole amplitude ji within a complex molecule. 
Let us assume that the value of /i is fixed so that ji gets effectively regularized only by the 
self-polarization cycles due to the presence of polarizable atoms within the host molecule, 
being a the total bare polarizability of the molecule. The net effect is that the spontaneous 
field emitted by the atom in the decay process gives rise to an induced dipole moment in the 
molecule which modifies the decay rate. If there existed other kinds of interactions between 
the fluorescent atom and the host particle, ji and k res might be modified by additional non- 
radiative effects. In the following, we will refer to the atom as emitter, to the molecule 
of polarizability «o as host particle and to the surrounding medium as host-medium. The 

2 — 

dipole moment of the system emitter- host-particle reads p — jl+ ^paQ u (f,r)fI, a being the 
renormalized polarizbility of the host particle. The perturbative series in Fig.8 allows us to 
build up the analytical expression for straight away from that of W w . The integration of 
the electrostatic part which yields the term proportional to a I must be removed with respect 
to the series in Fig.8 (a) for in the present case the source is the spontaneous emitter with 
transition amplitude fi. The rest of the terms remain unaltered but for the substitution of 
the induced non-radiative dipole moments €q(XqEq an d e o a oE* by the fixed dipole moments 
\x and jl* at the emission and reception sites respectively, 

r M - T H .(Mo) [ 1 _^3 ao + | , o2ao[27± + 7||] -H) 

-2u 2 |/i| 2 



3e c 2 ft |1 - ^kla + \kla [2 1± _ + 7 ,|]|2 



3{2 7 r- + 7||} ( 49 ) 



kn ^ r -.1 . k[ 



- ^{a }|-^ + 2 7± + 7|!| 2 J, (50) 
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We recognize again the power absorbed within the host particle in the last term and the 
power radiated into the medium in the remaining. Because any possible resonance of the 
host particle is assumed far from that of the emitter, we have ignored the regularization of 
Jt{2 T f}. 

(a) " V 

r ; im{o([Ot OWQ* £]♦ [Q + $#1 



FIG. 8: (a) Diagrammatic representation of Eq.(39). (b) Diagrammatic representation of 
Eqs.(49,50). 



IV. THE VIRTUAL CAVITY SCENARIO 
A. Virtual cavity vs. real cavity 

The 7-factors are the traces of the solution of the Lippmann-Schwinger stochastic equa- 
tion for the propagator of the self-polarization field of a point source, Q(f,r'), 

0ij(r, r) = G\f(f, - kl J d 3 rqGf (f, n) [(e km {n ~ r) - S mj ) |J Q mj {f^ (51) 
See [97] for an explanation about the notation. The brackets of the stochastic kernel 



efcm(/i — r) — S m j) denote average taking over all possible configurations of the sur- 

/ r,f - 

rounding host scatterers and both r and f are inside the emitter such that jr*— r\ < a <^ k^ 1 . 
In turn, the distinction between r and r 1 is just formal as it cannot be resolved. The explicit 
dependence of e(fi — r) on the emitter position f and the restriction symbolized by \^ 
signal the fact that f, r* are kept fixed inside the emitter as performing the average over all 
possible configurations of host scatterers. This way, any scattering process of the virtual 
photons in their way from and towards the emitter correlates to the emitter position. In 
case the emitter is in all equivalent to the rest of host scatterers within the host medium, 
such a correlation is the same as that among the host scatterers themselves. It is in this 
sense that translation invariance is only virtually broken and the exclusion volume around 
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the dipole emitter is referred to as virtual cavity. In case the emitter is distinguishable w.r.t. 
any other scatterer, that correlation differs. Hence, translation invariance is actually broken 
and the cavity is not virtual but a real cavity. The break of translation invariance, either 
virtually or actually, makes the computation of Q{r,r l ) different to that of the ordinary 
Dyson propagator. This is the reason why the propagator appearing in Eq.(22) for Fermi's 
Golden rule is not Dyson propagator as erroneously suggested in [35, 49]. In any case, it is 
important to emphasize that a translation invariant Dyson propagator does exist for normal 
modes in the virtual cavity scenario. In contrast, a strictly translation invariant propagator 
for normal modes does not exist in the real cavity scenario. 

This section and the next one are devoted to the computation of Q and the 7-factors in 
both the virtual cavity (VC) and the real cavity (RC) scenarios. Notice that not only the 
equivalence/inequivalence of the spatial correlation of the emitter matters in order to define 
each scenario but also the equivalence/inequivalence of the polarizability of the emitter w.r.t. 
that of the host scatterers. That is so because the 7-factors also depend on the polarizability 
of the emitter. In particular, the VC scenario requires that the dipole transition amplitude 
of the emitter satisfies the relation, «o = ^ifa-Io ■ «o being the electrostatic bare polarizability 
of the host particles. As pointed out in [50], the bare polarizability of a two-level excited 
emitter has the opposite sign with respect to that of the same particle in its ground state 
[44]. More importantly, the polarizability during the decay process is undefined. Therefore, 
stricto sensu only in the case of induced emission with all the scatterers of the medium in 
their ground state and for the case of a fixed dipole on top of a ground state host parti- 
cle, emission in a virtual cavity scenario holds. Hence, it is only in these scenarios that a 
translation-invariant actually polarized sourceless vacuum can be defined. 

The case of emission from a real cavity is much more generic and its treatment, at least 
formally, more complicated. There are nevertheless two situations where calculations sim- 
plify. These are, in the first case the excited emitter is weakly-polarizable and occupies an 
interstitial position in the host medium. In the second case, the interaction emitter-host 
medium is such that the emitter expels the surrounding host particles at a distance R much 
greater than the the typical distance between scatterers. 
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B. Computation of the virtual- cavity 7-factors 



Rather than solving the stochastic integral equation Eq.(51) directly, we will compute 
first the translation invariant bulk propagator G(r — r 1 ) which obey the Dyson equation, 

G l3 (f- f ) = Gf(f,r<) - klf d^Gf (r^rl) [(cfcm^i - r) - <^}] ^(rl, r<). (52) 

In contrast to Eq.(51), the average has no-constraint in this case. Hence, G(r — r f ) is a func- 
tion of r — r 1 for any two given points in the bulk. The field so propagated is the propagating 
coherent field or Dyson field, E D {r). That is, the field obtained by averaging point-wise the 
value of the classical Maxwell electric field E Max . (f) in the ensemble of spatial host scatterers 



configurations, E D (r) = \{E^ ax (f)}) -see Section VI Bl too. The photons which 

\ / ensemble 

propagate with G(f,r > ) are said on-shell or normal modes as they obey certain dispersion 
relations -see below. Because the set of fields {E^ ax (r)} are in one-to-one correspondence 
with the set of scatterer configurations and since each one obeys a linear Maxwell equation 

— * 

with certain dielectric function e m (r) specific of each particular configuration, Er>{r) obeys 
the same functional equation but with an averaged dielectric function ({e m (r)} 



ensemble 

Hence, Eq.(52) [11]. There exists a relation between the propagators G of Eq.(52) and Q of 
Eq.(51) which we proceed to find out. 

In Fourier space, isotropy allows to split the Dyson equation for G(k) in two uncoupled 
and mutually orthogonal algebraic equations, 

G ± (k) = G ( °\k) - k 2 Gf\k) x ±(k)G ± (k), (53) 
G l} (k) = G§>(k) - kl Gf{k) X \\{k) (?„(*), (54) 

where Gf\k) and G^\k) are given in Eq.(27) and x{k) is the electric susceptibility tensor 
such that e(k) = I + x(^) is the effective dielectric tensor. Statistical isotropy allows to 
decompose e and x i n longitudinal and transverse components, e± t \\(k), x±,\\(k)- For a 
medium made of a collection of disconnected dipoles, x±,\\ (k) can be expanded as a series of 
one-particle-irreducible (1PI) multiple-scattering terms of order n, 

00 

X Ul (k) = J2x ( :l(k)p n a n . (55) 

?1=1 

The functions x±,\\ are named 1PI as they contain integrals over undefined photon mo- 
menta which cannot be disentangled [11, 29]. In the above formula, p is the average 
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numerical volume density of scatterers and a is the renormalized polarizability of single 
scatterers. The functions X^(k) incorporate the spatial dispersion due to the spatial 
correlation within clusters of n scatterers. In particular, = = 1 incorporates 

all the self-correlation factors. In field theory terminology, x±,\\ are proportional to the 
photon self-energy functions, H± t \\(k) = — &oX_i_,||(A;) -not to be confused with the dipole 
self-energy. Alternatively, Eqs. (53,54) can be written in terms of the so-called T-matrix, 
Tyi(fc) = E ±> ||(fc) + Z U] (k)G U] (k)X U] (k), as 

G U] (k) = G%(k) + G%(k) T U] (k) Gf\(k). (56) 

Two remarks are in order at this point. The first one concerns the topology of the host 
medium. Because the emitter is treated as a point emitter and it is equivalent to the host 
scatterers, a cermet topology is inherent to the virtual cavity scenario. In the expression of 
Eq.(55), spatial dispersion in x cannot be disregarded as there must be at least a minimum 
distance £ between scatterers which determines their exclusion volume. Hence, this distance 
is the virtual cavity radius. For the sake of consistency with the linear medium approxi- 
mation, £ must be greater than the dipole radius to avoid overlapping between orbitals and 
so that the dipole approximation be valid. Second, in view of Eqs. (53,54), longitudinal and 
transverse modes of the Dyson field do not couple to each other as traveling throughout a 
random medium the same as it occurs with photons in free space. Eqs. (53,54) can be solved 
independently yielding the Dyson propagator components, 

1 



G ± (k) 



G ^ k) = wrmv (57) 

In terms of free propagators and self-energy functions they can be depicted perturbatively as 
in Fig.9. A more detailed examen in Section VI will show that longitudinal and transverse 
bare photons -i.e., the ones with propagator G^- do couple necessarily when they experience 
multiple scattering processes. In Eqs. (53,54), longitudinal and transverse bare photons enter 
both x±(k) an d Xll(^) by means of the spatial correlations among scatterers. 

Getting back to the computation of Q and the 7-factors, we first observe that 27j_ and 7y 
are made of 1PI diagrams in which the starting and ending points coincide at the emitter 
location. Those diagrams belong to x(h) and amount to the so-called recurrent scattering. 
As a matter of fact, the emitter itself enters the bulk propagator of Eq.(57) as an ordinary 
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FIG. 9: Diagrammatic representation of the Dyson propagator C 

scatterer in the VC scenario. The second observation is that, as the averaging process in 
Eq.(51) is subject to the fixed location of the emitter, every scattering event in the 1PI 
diagrams contributing to 7 is correlated to the emitter either at the near end or at the 
rare end of each diagram, indistinguishably. This is a consequence of reciprocity. Taking 
advantage of this feature in every multiple-scattering diagram like that in Fig. 10 (b), we 
can attribute all the irreducible correlations of the intermediate scattering events to the 
emitter on the left. By so proceeding, we end up with an effective separation of all those 
pieces irreducibly correlated to the emitter on the left completely disentangled from those 
non irreducibly correlated which form non-lPI pieces on the right. The sum of the 1PI 
pieces on the left amounts to x/p<3, where the factor 1/ pa stands for amputation of the the 
first random scatter which enters the diagrams of x 111 favor of the 'virtually' fixed emitter 
location. The sum of the non-lPI pieces on the right amounts to the bulk propagator G. 
Therefore, we end up with the formulae [48] -Figs.l0(c, d), 

g vo (k) = _L x±{k) a±{k) = _L gp^gpp , (58, 

grik) ^± &xt(k)Glm ^_^_. (59) 

The above expressions contain both G^ and G^ which carry divergences. The expressions 

vc _ r d*k^ r 2 X ±(k)/(pa) _ 2G m (k) ] m 

7± " J (2tt)3 L*g[l + Xj .(Jfe)] - ife» 2G±{k) \> m 

7 " -J (27r)3 [pa kl[l + *„(*)] G H {k) \> (61) 
are however fully convergent. We emphasize that the above expressions for the virtual- 
cavity 7-factors are exact under the assumption that the emitter can be treated as a point 
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FIG. 10: (a) Feynman rules. Only two-point irreducible correlation functions have been used 
for the sake of simplicity. (6) Diagrammatic representation of the equivalence between multiple- 
scattering processes amounting to Q. (c),(d) Diagrammatic representations of Eq.(60) and Eq.(61) 
respectively. 

dipole equivalent in all to the rest of host scatterers. The explicit inclusion of a -as given 
by Eq.(42)- in the computation of through Eqs. (58,59) makes it depend implicitly on 
their integral quantities 7_q| given in Eqs. (60,61). Later on we will see that the form of the 
above expressions for 2^°' and ^ c will allow us to classify the nature of the emission in a 
transparent manner. 

Alternatively, we can write Q vc (k) in other forms making use of Dyson's equation. In 
function of the bulk propagator of Eqs. (53,54) it reads 



k^pa 



G 



(o) 



In function of the T-matrix of Eq.(56) it reads 



-i 

k\pa 



T ul (k)G%(k). 



(62) 



(63) 



Finally, we write the Lippmann-Schwinger equation of Eq.(51) for the self-polarization 
local field propagator in Fourier space, 

Ql c {k) = Gf\k) + Gf(k) E v ± c (k) gl c (k), (64) 
QY c (k) = Gf\k) + G<°> (*) Ej c (k) g™( k ), (65) 
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where 



Sift*) = -^)77T [' '^P • *&a, B (*) G^Wj ( 66 ) 

X-L,||(*) _L,||( ) P " J 

is the stochastic kernel expressed as a function of x, and pa. 

Because the transference matrix (i-matrix in brief) formalism is profusely used in nu- 
merical simulations which involve discrete configurations of point dipoles (see eg. [51, 52]), 
we will next show that analogous formulae to the ones above are obtained for a fixed con- 
figuration of point dipoles. Let us take the m th configuration of the statistical ensemble 
described in Section II C. Its corresponding self-polarization vacuum state there was de- 
noted by \4> m ) s ' p '- Let us assume the configuration consists of N + 1 host scatterers with 
fixed position vectors denoted by {f;} with % = 0, N. We then decide to excite the dipole 
at r*o with an external monochromatic field of frequency u = k$c such that E exc (f) = Eft 
iff r — fo and E exc (f) = otherwise, as in Section III A. All the dipoles are equivalent and 
have renormalized polarizability a. Therefore, the dipole moment of the emitter reads 

p(f ) = e aE%. (67) 

The total field at the emitter location is the sum of the exciting field plus the self-polarization 
field which the dipole moment p(f ) itself creates at r . That is, 

E(f ) = E^ + ^g(f ,f )-p(f ) 

= E% + k 2 ag(r ,r )-E%. (68) 

The second term on the r.h.s. of Eq.(68) is the field which propagates freely in space from 
all the induced dipoles to r , 

k 2 ag(f , f ) • Eg = M V G(°)(f , n) ■ m- (69) 

Except for the emitter, all the rest dipoles are only induced by their mutual interactions. 
The t-matrix, with components t(fi,fj), in the same spirit as that in Eq.(56), yields the 
dipole induced at some point fj as a result of its interaction with the dipole excited by an 
external source at some other point fj. That is, 



X-L,||(*0 



2., /u\rt°) 



m = (f)' 1 ?i) ■ ZeUrj). (70) 



{rj} 
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Analogously to the susceptibility function in Eq.(55), t^(fi,fj) can be expanded in powers 
of — k^a, t(fi,fj) = J2 n= ii^ n \ri,fj). Each term in the sum contains n factors — k^a and 
n — 1 tensors &°\r m , f{) which propagate the field though all possible paths connecting the 
points r*j and fj. As an example, 

= (~*fc) 4 E ^(n,^) • G (0) (w,) ■ G< >(*5,3), (71) 
{n,f m } 

where restrictions apply over the indices of the sums to incorporate correlations. The analog 
to the T(k) matrix of Eq.(55) is just 

f(k,k') = E tfafye'fc*-*^. (72) 

{n,rj} 

Because in our case the only externally excited dipole is that at r , Pi = — tffi, tq)-Eq. 

Inserting this formula into Eq.(69) we obtain 

Q(fo, r ) = — ^ E ° (0) ^o, n) ■ t(n, ro), (73) 
which is the analogous expression to that in Eq.(63) but for a fixed configuration of scatterers. 



V. THE REAL-CAVITY SCENARIO 



The computation of the 7-factors in the real cavity scenario is generally quite more com- 
plicated as an effective medium cannot be strictly speaking defined and translation invariant 
bulk propagator and suceptibility tensors do not exist. There are two main differences with 
respect to the virtual cavity scenario. The first one is that the emitter does not behave 
as an ordinary scatterer of the host medium regarding photon propagation. The second is 
that the presence of the cavity induces additional correlations among the scattering events 
experienced by photons near the cavity surface. That is, because scattering events are corre- 
lated to the emitter through the stochastic kernel of Eq.(51), they indirectly correlate among 
themselves beside the inherent spatial correlation of scatterers in absence of cavity. 

There are nevertheless two situations in which approximate solutions can be found. The 
first one corresponds to the case in which the emitter hardly alters the host medium. For this 
to be the case, the polarizability of the emitter must be much weaker than that of the host 
scatterers and its cavity radius R much less than the typical distance between scatterers, 
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R <C p -1 / 3 . In his situation the foreign emitter is said weakly polarizable and interstitial. 
The second situation which can be treated analytically is that in which the cavity is large, 
R S> p -1 / 3 . The density function of host scatterers is altered in a large patch but remains 
unaltered locally within distances of the order of £ where correlations matter. In this case, 
the foreign emitter is said to be substantial. 



A. Interstitial weakly-polarizable impurity within a small cavity 

The two-point correlation of the host scatterers to the foreign emitter reflects in this 
case on a slight modification of their average density function in the neighborhood of the 
impurity. The radial distribution density function with respect to the emitter site reads 

p{r i )=p[l-Q(r i -R)\, R « p~ 1/3 , (74) 

where r$ is the position vector of a generic scatterer with origin at the emitter location 
and p is the 'would-be' uniform numerical density of host scatterers in absence of emitter. 
As an example, we depict in Fig. 11 (a) the diagrammatic representation of the scattering 
amplitude for the process in which a photon emitted by the foreign emitter is scattered by 
a host scatterer and absorbed by another one, J-'^(k), 

pV){k) = p(pa) 2 j d 3 nd 3 r 2 e^ 2 [Tr{G (0) (ri) • G (0) (r 2 - n)} 

x [l-e(r 1 - J R)][l-e(r 2 - J R)][l-/i(|f 2 -f 1 |)], (75) 

where h{\r 2 — T\\) is the two-point correlation function between host scatterers with support 
within a correlation volume ~ £ 3 . It is clear from Fig. 11 (a) that in general it is not possible 
to disentangle the correlation of the emitter to each scatterer event. However, for the case 
that £; £ ^ 1; k R <C 1, one can apply the overlap approximation [15] for those modes 
k <C This is so because those modes cannot distinguish between the exclusion 

volume of the host scatterers and the real cavity. That is what Fig. 11 (b) depicts. In principle, 
the overlap approximation in this scenario demands that the only relevant correlation in the 
medium be given by the Heviside function of exclusion volume. This is clearly the case 
of a weakly correlated host medium (eg. a diluted gas like that of Fig. 12(a)) for which 
h(\f 2 — fi\) ~ — 0(|r 2 — f*i| — £) with £ <C p -1 / 3 . For a highly correlated medium with 
£ ~ p -1 / 3 (eg. a solid like that of Fig. 12 (6)) a term proportional to 5^(r — £) amounting 
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for first neighbors must be included in h(\f 2 — fi|). Nevertheless, for /c £ <C 1 that term is 
irrelevant in comparison to the exclusion volume function and the overlap approximation 
for k <C is equally valid. Therefore, the condition k t; <C 1 is generally sufficient for 
application of the overlap approximation. In addition, this situation corresponds to that 
in which an effective medium exists with Xeff the constant effective susceptibility valid for 
those /c-modes with h£ <C 1. The contribution of transverse modes with k£ <C 1 to LDOS, 
I\ and W± is well approximated by the formulae of the VC scenario which we will compute 
in Section VI B 3. The corresponding 7-factor will be found to be 

^^^SS?ii G ' l ' (k) - (76) 

For those modes k > the equivalence with the virtual cavity formulae Eqs. (60,61) is 
only possible for the case that -R ~ £ -see Fig. 11(c). This is only applicable to a weakly- 
correlated host medium as otherwise R would not satisfy the interstitial condition R <C p -1 . 
That is the case of the Maxwell-Garnett dielectric which will be studied in Section VII A. 
This has the nice implication that, if the impurity is an excited host scatterer of a diluted gas, 
the r y vc factors computed in Section IV B are applicable in good approximation provided 
that the polarizability of the excited atom during its decay process is weak in comparison 
to that of the rest of scatterers. 

B. Substantial impurity within a large cavity 

As mentioned above, the presence of an emitter which substitutes one or several host 
scatterers can modify sensitively the homogeneity of the host medium. In particular, this 
is the case of a host scatterer which is promoted to an excited atomic state in a strongly 
correlated medium. In comparison to the scenario of the previous Subsection, this is the 
same situation but with £ ~ p -1 / 3 . Because the polarizability of the excited scatterer is 
different to that of the rest, namely a exc- , it behaves as a foreign emitter which occupies an 
effective volume of the order of p -1 . Thus, the approximate average polarizability density at 
a point r away from the emitter is given by pa[l — @(r — p -1 ^ 3 )] + pa exc Q(r — p -1 / 3 ). This 
makes the electrical susceptibility to be clearly not a homogeneous function. For instance, 
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FIG. 11: (a) Diagrammatic representation of Eq.(75). (6) Overlap approximation restricted to 
those modes fcC^ 1 subject to fco£ <Cl. (c) Overlap approximation for shorter wavelengths valid 
only for i? £. 

the two-scattering term would read 

x "a[l - 0(r 2 - p- 1 / 3 )] + pa MC .0(r 2 - p- 1 / 3 )] - r^MI^ - ^0.(77) 

That implies that the susceptibility function evaluated at r 1? r 2 close to the emitter can be 
very different to that at fi, r 2 far from it. 

However, whatever the nature of the foreign emitter is, if the cavity radius is large in 
comparison to the mean distance between scatterers, R ^> p -1 ^ 3 , the medium will look 
statistically homogeneous at scales of the order of p" 1 / 3 as seen from the emitter. This 
is so because the electric field which propagates from the emitter is statistically uniform 
within the interval [R, R + p -1 ^ 3 ] for R ^> p -1 / 3 . Thus, the average value of the numerical 
density of host scatterers at a point r*j away from the emitter reads in good approximation, 
p(fi) ~ p[l — 0(rj — R)\. Because the typical correlation length between host scatterers is 
£ < p -1 / 3 <C R, the cavity exclusion volume factor — 0(r — R) is common to all the scatterers 
which are correlated among themselves and enter the susceptibility function. This has the 
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nice implication that solids with well-defined susceptibility function x(k) can De treated 
this way with no-need to know the individual polarizabilities of their constituents. Any 
fluctuation not yet considered for the case of point-host-dipoles (eg. excitons) is assumed 
integrated out in provided that the wavelength of those fluctuations is much shorter 

than R. In the following we develop a formalism which bases on the definition of quantities 
analogous to x(k) an d G{k) but for a non-translation invariant medium. 

Let us formulate mathematically the above approximation. Let g(r) be the full two-point 




FIG. 12: Two dimensional sketch of the real cavity scenarios which can be treated analytically in 
some approximation, (a) Weakly correlated host medium (eg. gas) with an interstitial impurity 
[x. (b) Strongly correlated host medium (eg. glass) with an interstitial impurity, (c) Large cavity 
enclosing a substantial impurity. 

correlation function of the emitter with the host scatterers, 

g(r) = l + h c (r), (78) 

where hcij-) is the irreducible piece. In general, because the cavity is macroscopical, it can 
distort the homogenous distribution of host scatterers in the surrounding medium. In such 
a case hcij-) must contain a maximum at r ~ £ that makes account of the overdensity of 
first neighbors. For the sake of simplicity, we will restrict he to a step function which stands 
for the cavity exclusion volume, h c (r) — 0(r — R). The Fourier transform of g(r) reads 

g(k) = (2n) 3 5 {3 \k) + h c (k). (79) 

Because hc(r) has support in r < R, it is therefore expected that its Fourier transform gets 
support in k < 1/R. On the other hand, in case that spatial dispersion in e(k) be relevant, 
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such a dispersion must be of the order of l/£. Therefore, any convolution of g(\k — k f \) with 
the self-energy operator £(&/) and any generic function f(k') can be approximated by 



/ 



^ - k\Mk) ■ f(k) » E(fc') • | ^<KI# - *!)/(*)■ (80) 



(2tt)3 ; 

In Eq.(80), the difference between k and fc' is negligible in comparison to l/£ in the range 
of momenta where hc(\k — k'\) takes nearly constant value and S(/c) is not zero. In turn, 
that implies that the correlation of the emitter with any 1PI (multiple)scattering process 
can be approximated by the correlation of the emitter to any of the scatterers involved in 
such a process. In particular, things get mathematically simpler if the correlation functions 
connect the impurity either to the first or to the last scatterers within the 1PI diagrams of 
E(/c). This allows to write the series for the RC propagator, Q RC , as two apparently different 
expansions, (a) and (b), as depicted in Figs.l3(a, b) respectively. The corresponding series 
read 

? w _ c}" "' g lir t~"-' 

n=0 n=0 

oo oo 

g™ = gf + £ ^ (n) , Qf c = $T + E g T ] ■ ( 82 ) 

n=l n=l 

In each series, the n th terms read respectively, 



: a ~ ! ■+ \ ; + \£ .., \ ; + ;/ \ "\ ;+ 
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(b) 



FIG. 13: Diagrammatic representation of the self-polarization propagator Q according to 
Eqs. (81,83,84) -series (a)- and Eqs. (82,85,86)) -series (&)- in the large cavity scenario. 

g a ± {n \k) = Ki\k)E ± (k)Gf{k), (83) 
g< n \k) = tf\k)^(k)G§\k), (84) 
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and 



g b l n \k) = K { r 1 \k)i; ± (k)Gf\k)i; ± (k)K ( l\k), 



(85) 
(86) 



The recurrent formulae for the partial pseudo-susceptibilities -the reason for this nomen- 
clature will get clear later on- read 

d 3 *;' 



g(\k'-k\) (l + cos 2 fl) 



K^\k) 



(2tt) 3 - 

x (fc')Ex {k')G < f > (k 1 ) 

+ sin 2 6 ^(k'^ik'jG^ik') 
d 3 k' 



for n > 1, 



(87) 



-MW-k\) 



cos' 9 



(2tt) 3; 

x ^(k'^k^ik') 
+ sin 2 9 K±~ 1 \k')E±(k')G ( ± (k') for n > 1, 



(88) 
(89) 



and «^||(ife) = -PE^i^A;)- 1 = l/x±,\\(k) for n = 0. 

In Eqs. (87,88), Ej^^/c') can be factored out of the integrals as E_q|(A;) in application of the 
approximation in Eq.(80). Therefore, the building blocks of the series are G { ^(k) and the 
irreducible pieces of K^(k). Hereafter we will refer to the 1PI parts of K^(k) as cavity 
factors, 



C±(k) = - I d 3 r e*\(r)Tr{G (0) (r)[I -k®k}} 
-h c (\k'-k\)\Gf\k') 



2 _ 
1 f d 3 k' 



2 J (2vr) 3 
+ (A;')cos 2 ^ + G\°\k') sin 2 9\, 

C\\(k) = J d 3 r e ilr h c (r)Tr{&°\r)[k <g> k}} 

r a^u' 

' ;h c (\k'-k\) 
J l Z7 ' ) 
f „(n\ , , , 



(90) 



' (2tt) 3 

Gf^^cos 2 ^ + (fc')sin 2 0], 



(91) 



and so «^„(fc) = -k 2 [G { l\(k) + C ± ,,|(A;)]. 

Also, Eqs. (85,86) resemble the formulae of the 7-factors for the virtual cavity scenario 
-Eqs. (58,59). That is, if we define the total pseudo-susceptibility as K±^\(k) = X] n =o K ± \(k) 
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we can write 

Gf C {k) = n ± (k) X ±(k)G^\k), (92) 
gf c {k) = Kl] (k) Xl] (k)G$\k). (93) 

However, the above equations are not yet quite similar to those for the VC scenario. In 
particular, the propagator on the r.h.s. of Eqs. (92,93) is that in free-space whereas it is 
the bulk propagator in the virtual cavity formulae. The reason being that K± i \\(k) contain 
both 1PI and non-lPI processes. We can separate n±,\\(k) into 1PI and non-lPI (N1PI) 
pieces, Kj_ t \\(k) = K^(k) + k 1 ^ 1 ^ 1 {k) according to the following decomposition in partial 
pseudo-susceptibility functions, 

^l]\ I{n \k) = 4"- 1) (A;)S ±i ||(A;)G'f (A;), n > 1, 
^ n \k) = _ K™ w <»>(]t), „ > 1, 

*^(k) = K%(h), <r (0) W=0. (94) 

Using this decomposition, we can write Eqs. (92,93) in the form, 

Qf C {k) = ^ ± PI (k)x±(k)G ± (k), (95) 
Gf c {k) = «jj w (fc) X ||(fc)G||(fc), (96) 

where G± i \\(k) is the Dyson propagator given by Eqs. (53,54) for a 'would-be' homogeneous 
medium in absence of the real cavity. As a matter of fact, in passing from Eqs. (92,93) to 
Eqs. (95,96) one can use the same arguments as those employed in the virtual cavity scenario 
to push the effective bulk propagator to the right as in the diagrams of Fig. 10(6). Thus, 
the above equations resemble those expressions for the VC scenario in Eqs. (58,59) with the 
replacement of by (k)x±,\\(k) . Therefore, the relation between the propagators 

G±^\(k) and G±9\(k) for two identical random media which differ just by the presence/absence 
of a real cavity enclosing the emitter is given by 

G RC (k) 

0^ = pa«5f(*). (97) 
Finally, we write the Lippmann-Schwinger equations for the large real cavity scenario, 

G±°(k) = G ( °\k) + Gf ^Zfc^g^ik), (98) 
gf c {k) = Gj°)(fc) + G^\k) Ef(k) Gf c {k\ (99) 
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where the stochastic kernel is 



Sfji(Jfe) = ^ 1 - K^Uk) X ± \\{k) + k 2 oX ± \\{k) G { l\(k) . (100) 

VI. THE NATURE OF DIPOLE EMISSION AND THE COHERENT VACUUM 

In this section we first study a simple model in which dipole emission can be studied 
microscopically in full detail. We compare the results with the usual Lorentz-Lorenz (LL) 
and Onsager-Bottcher (OB) formulae and other more recent works. Next, the same decom- 
position is performed over the general formulae of the previous sections. Our microscopical 
approach allows us to distinguish between the transverse and the longitudinal, the coher- 
ent and the incoherent, the direct-coherent and the induced-coherent, and the dispersed- 
incoherent and the absorptive components of the dipole emission. Special attention is paid 
to the phenomenon of radiative/non-radiative energy transfer. Our study reveals why the 
nature of dipole emission was erroneously interpreted in previous works. 

In connection with QFT formalism, we define a coherent vacuum associated to coherent 
emission. We find its relation to the sourceless vacuum in a random medium, which is itself 
associated to direct emission. 



A. Decomposition of dipole emission in the single scattering approximation 

Our simplified model consists of a stimulated emitter embedded in a diluted host medium 
made of spherical inclusions of electrostatic bare polarizability a®. The scatterers present a 
minimum spherical exclusion volume of radius £ and density p <C £~ 3 . The emitter is excited 

— * 

by an external field Eq whose frequency, u = cko is such that k £ < 1. A similar model 
was studied in [23, 53, 54]. The single scattering approximation implies that virtual photons 
experience, at the most, a single scattering event with the scatterers in the host medium. 
This way, the only relevant correlation is that given by the two-point function h(r) = — 0(r — 
R). If the cavity radius R satisfies R = £ and the electrostatic polarizability of the emitter 
is the same as that of the host scatters, ao, this set-up corresponds to the VC scenario. It 
corresponds to the RC scenario otherwise. Let us assume without loss of generality that 
R = £ and let us take a e as the electrostatic polarizability of the emitter which, eventually, 
can be taken equal to a . This implies that, in the single scattering approximation both 
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FIG. 14: (a) Feynman's rules . (6) Diagrammatic representation of W^. The first three diagrams 
stand for those terms in Eq.(103). The next three diagrams take account of the self-polarization 
terms of Eqs.(119,120). 

scenarios are equivalent with just exchanging the cavity factors appearing in the RC scenario 
of Section VB with the two-scattering term of the susceptibility in the VC scenario, yf® ■ 
They relate through x±\ — —ko(pa) 2 C±^\ R= ^, where a is given by Eq.(42) and C±^(k) are 
given in Eqs.(90,91). The susceptibility of this model reads xfj.W = P<5 + X|| 2 j_(^)- I 11 the 
following, the superscript [1] will denote single-scattering approximation. 
The quantity to compute is 

W " = tf'V... Jw'l 9 I m 9 ^" + I 2 ' < 101 > 

oc |i + _A;^ e [2 7 ^ + 27 ± + 7||J| 2 
Pictorially, 7M is given by the second and third diagrams of Fig. 14(6). It reads, 

7 [1] = Tr{ J d 3 n G^i^i-klp^G^in) [1 - e(n - £)]}• (102) 

Let us disregard for now self-polarization cycles on the emitter. In doing so, the emitter 
cavity is said empty. In the limit fco£ «C 1, up to O(0) terms in A; £ we get 

w W,em P ^o| 1+ <^ {p ~ } + + ^] }, (103) 
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where W° = ^^r\E£ | 2 is the power emitted in free space. The 0(0) terms of Eq.(103) 
equal the usual Lorentz-Lorenz (LL) and Onsager-Bottcher (OB) formulae in absence of 
back-reaction effects [55-57] with = 1 + pa, 



W, 



[i] 



LL 



w, 



[l],emp 



OB 



= W»(l 



7 



2K{ e [ 1 l} + l 



-K{pa}). 



9fc{VeW}, (104) 



(105) 



Same results as those in Eqs. (102-105) were obtained by the authors in [23] using a Quantum 
Optics formalism. The factors ^ M^±g j anc i are referred to in the literature 

as virtual cavity and empty cavity local field factors respectively. The factor A/K{e} is 
sometimes refereed to as bulk factor. Both Wll and W^f are commonly attributed in the 
literature to transverse emission. Our next task is to clarify the actual nature of the terms 
in Eqs. (102,103) through a rigorous microscopical analysis. In Fourier space, Eq.(102) can 
be decomposed into transverse and longitudinal parts, 

d 3 k 

~G { l\k')5^(k' -k)Gf\k) 



2 7 1_ 1] = Ipa 



Gf\k) + C ± (k) G K l\k) = 2pa 



,(0)/ 



& A k 



d 3 k' 



(2tt) 3 J (2tt) 



+ -Gf\k')h(\k' - k\)(l + cos 2 e)G^>(k) 



(o), 



- ^Gf{k')h{\k> - k\) sin 2 6G^\k) 



(106) 
(107) 
(108) 



pa 



tfk 



G\°\k)+C\\(k) G K l\k)=pa 



i(0), 



x 



(2tt) 3 

Gj 0) (ifc , )* (3) (fc , -fc)G f [| 0) (A;) 
+ G^AOMI*' " fc|)cos 2 0Gjj O) (ife) 
+ Gf\k')h{\k> - k\) sin 2 Cf 0) (/fc)" 



where h(q) = [sin (g£) — g£ cos (g£)] and we can write, 



d^A; 



< ] = W2[1-^{2 7 W}], 



[i] 



2ttc 



d 3 A;' 



(27T) 3 ,/ (2tt) 



(109) 
(110) 
(111) 



112) 



In the integrals of Eqs. (106-111), the only imaginary pieces are those associated to the 
poles of G < f ) {k) factors. Therefore, ^{tIui} is composed of terms of the form $t{pa} times 
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residues computed at k = k plus terms of the form Sjpo;} times the real part of the 
integrals. Alternatively, we can write 27^' and 7m 1 ' in terms of spatial integrals, 



,3 r 



2 7 i 1] = p«Tr{ J d 3 r[G^.(r) + G^ at (r)] ■ G^.(r) [1 - 0(r -£)]}= p& J d 3 

x K{G { l.(r)-G^.(r) (113) 

" GS.(r)-Gi(r)6(r-0 (114) 

- Gi°i,(r)-Gil(r)0(r-O}, (115) 

7™ = p«Tr{ J d 3 r[G[t.(r) + Gfi,(r)] ■ G^,(r) [1 - 9(r - $]} = pa J d 3 r 

x Tr{Gil(r)-Gil(r) (116) 

- G2i,(r)-GSi,(r)0(r-O (117) 

-0}. (118) 



- ^l(r).GS,(r)9(r 



The correspondence between the integrals in Fourier space and those in spatial space is 
obvious. However, such an obvious correspondence is only possible in the single scattering 
approximation. This is so because the bulk propagator is identical to that in free space 
but for the substitution k$ — > k$ so that it does not present spatial dispersion and 
there exists a direct identification between the radiative and electrostatic bulk propagators 
in spatial space and transverse and longitudinal propagators in Fourier space respectively. 
Note also that, in absence of the correlation function, there is no coupling between radiative 
(transverse) and electrostatic (longitudinal) modes. 

Next, we search for the terms in Eqs. (106-118) which amount to those of Eq.(105). The 
free space term 27^ together with that in Eq.(106) (or Eq.(113)) yield W°(l + |3?{po;}), 
which corresponds to the bulk factor W^°^{VeW}. This is the emission directly radiated by 
the dipole into the medium -see Fig. 16(a). Eq.(106) (or Eq.(113)) corresponds also to the 
diagram of Fig.5(d 2 ), which explains why of the two transverse propagators there. Eq.(107) 
(or Eq.(114)) yields higher order terms. The common terms to 27^' and 7m 1 ' are those in 
Eq.(108) (or Eq.(115)) and Eq.(lll) (or Eq.(118)) which amount to W^{pa} each. They 
equal one local field factor each, ^ Mllilij j or ^_^||iL.^_ The first thing we learn from 
this analysis is that, contrarily to the common assumption, one of the local field factors is 
associated to longitudinal emission while the other one together with the bulk factor belong 
to transverse emission. The latter is associated to the poles of the bulk transverse propagator 
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of the general expression in Eq.(58) and form part of the coherent radiative emission -see 
next subsection. On the contrary, the former is associated to the imaginary part of x\\ which 
contains a transverse bare propagator (see Eq.(91)). It is part of the incoherent emission 
which is dispersed. The existence of the common term in Eqs.(108,lll) is a consequence 
of reciprocity. However, while it gives rise to coherent emission when the coupling reads 
G\\ — G±, it gives rise to dispersion as read in opposite direction as in Fig. 15(6). 

Regarding the ^-dependent terms of Eq.(103), they both are proportional to ^s{pa}. The 



\ h(r) = -0K) h . 

-= = C * ^ r 
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M' 
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FIG. 15: (a) Feynman's rules, (b) Common vertex derived from Eqs. (108,111) which couples 
longitudinal to transverse modes, with \k x k'\ = sin (9. (c) Longitudinal-longitudinal vertex from 
Eq.(llO), with \k • k\ = cos8. (d) Diagrammatic representation of the origin of the coupling of 
single scattering radiative corrections to longitudinal external modes. 



leading order one, is given by Eq.(109) plus Eq.(llO) (or Eq.(116) plus Eq.(117)) and 

is associated to the imaginary parts of both the longitudinal bulk propagator -Eq.(109)- and 
the longitudinal susceptibility -Eq.(llO). The term has two identical contributions 

of value | coming from the common term to and ^ in Eq.(108) and Eq.(lll) 
respectively. They are associated to the imaginary parts of x±(k) an d X\\(k) respectively. 
All these terms proportional to S{o} are related to absorbtion in the host scatterers. 

Finally we consider the back-reaction of the host medium on the emitter with electrostatic 
polarizability a e . To this aim and at first order in p we have to compute the self-polarization 
diagrams in the second row of Fig. 14(6). For the sake of simplicity, in the following we 
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disregard absorbtion and radiative corrections in a and keep leading order terms in a and 
k £. That is, ~ 1 + pa is fully real and 

^iM + I^^i^fi- (119) 

Plugging Eq.(119) into Eq.(lOl), we obtain 

W^^wg ] r P [l + ~pao], (120) 

where Vg is the volume of the cavity. One can verify the agreement of the above expression 
with the OB formula which incorporate self-polarization terms at leading order in ao, a e 
[14, 57], W^"" = ^BRPJ^fW*)'*^- 

B. General decomposition of dipole emission 

Next, we go beyond the single scattering approximation. Again, let us consider Eq.(lOl) 
in absence of absorbtion in the emitter. In order to examine the nature of the emission, it 
is convenient to write the 7-factors in terms of the propagators of both the virtual and the 
real cavity scenarios using Eqs. (58-59) and Eqs. (95-96) respectively, 

j3, 



121' 



u/vc f d k X±{k) ^ cm-X||(*0 n (U\\ 

W ™ K / ^[ 9 { 2 «^V(*)G ± (fc)} + 9f{«} w x||(*)G H (*)}]- (122) 

From now on, we will just work with W^ c bearing in mind that the computations are in 
all equivalent to those for but for the replacement x±,\\ ~^ P&k 1 j P \\X±,\\- F° r the sake of 
simplicity we will drop the scripts VC and uo and we will denote the proportionality constant 
omitted in the above equations by W Q . 

The first obvious decomposition is that between transverse, 2W±, and longitudinal emis- 
sion, W\\, 

2 ^ = H, °/(05 29{ ^ GlW> ' (123) 
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1. Coherent and incoherent emission; direct and induced coherent emission 



Coherent emission is that associated to the coherent propagating field. That is, the 
emission whose modes satisfy the same dispersion relations as the normal modes which 
propagate through the bulk [12, 36, 37], 

k%e±(k) — k 2 \ k=k p™p- = for transverse modes (125) 

and €\\(k)\ k=k pro P . = for longitudinal modes. (126) 

These are the poles of G±(k) and G\\ in the integrands of Eqs. (127,129). Therefore, we 
identify from Eqs.(121,122), 

^- 3 2^{^l}Q{G ± (k)}, (127) 
{Ztt) 6 pa 

29{M?W-l(*)}, (128) 
{Ztt) 6 pa 

n xm mG]l{k)h (129) 



2 \yCoh. 


= w 


2W[ ncoh - 


= W 


w Coh. 


= W 


yylncoh. 


= Wo 



(2tt) 3 pa 

^ { m mGm . (13 o) 

Because longitudinal normal modes need of material support to propagate and are subject 
to very specific arrangements of scatterers [42], W^ oh ' can be ignored for practical purposes. 
Let us examine firstly W Coh " in the framework of Classical Optics. It is given by 

3 r 

W Coh - = ^-JPo\ 2 J dVTr_{$ft{x(f-0/P«}-^{G'(^r;c i ;)}}, (131) 



where both x an d G are written in the spatial space representation for convenience and 
p = eoaE^{r) is the dipole moment induced by the external field E^{r) on the emitter. 
Considering the fields classically, the fluctuation-dissipation relation reads 

9{G w (f ,0} = -^|(^(0 K*(r)), (132) 

where E^if) is the u;-mode of the coherent-Dyson field and the script D stands both 
for Dyson and for direct emission. Using Eq.(132) and writing 9?{x(r — r')/pa} = 
[&{x(r- ?)/ pa} - 5^(f- r*)I] + 6®(r - F)I in Eq.(131), we separate explicitly the field 
emitted directly by the source dipole from that which is emitted by the induced surrounding 
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Qh 

7TUJ , 

+ 6h< 



dipoles, 

w Coh. = - lfQl 2 {E . { ^ . ££* (r -^ (133) 

■|p | 2 Tr_{ J dV^[K{x(f-^)/pa}-5( 3 ^f-^)I] •^(r y )^*(f)^}(134) 

The term in Eq.(133) is the coherent power carried by the field directly emitted by the 
source dipole into the bulk as if it were an external source. In Eq.(134) we can identify the 
field emitted at r* by the induced dipoles sited around the source which propagates towards 
the source located at r, 

E*{f) = J dV[9fc{x(r - r )/pa} - 5 (3) (f - f*)l] ■ E%{?), (135) 

where the subscript / stands for either induced or indirect. Therefore, we can write 

w cah. = w cah. + W cah. = ™^2 ^\E^(f)\ 2 } + $t{(E?(r)-E<g(r))}]. (136) 

As expected, the first term is the coherent power emitted directly by the dipole source while 
the second term is the coherent power associated to the interference of the field emitted by 
the source and that emitted by the induced dipoles. 

As an example, let us consider wf^' Coh " in the single-scattering model of Section VIA. 
It is given by Eq.(115) -modulo the appropriate prefactors- and its diagram is that of 
Fig.l6(6). Making the identification [3?{x [1] (f - f*)/pa} - 5 (3) (f- r*)I] ~ k%pa Q(\r - r*| - 
^KjG^ (r — r 1 ')}, and applying the fluctuation-dissipation relation, it reads 

Wl 1] ' c ° h - = C^Tr{ J dVGSL^-O ' im?-?\-Oeopao(EZ(7)}E^)}, (137) 

where the Dyson field is fully transverse. The quantity within square brackets is the 
polarization density induced by the coherent field on the surrounding host scatterers, 
p°(r r ) = 0(|r — f| — ^eopaoE^f 1 ). The electrostatic components of the fields emitted 
by the induced dipoles within the sphere of radius £ add up at the emitter location, r, in 
the form of and induced effective field, 

E'nr) = ^ / dVG2L.(r - f) • (138) 



It is the transverse component of the above field that interferes with E D (r), 

w m, Coh . = w o^ { g r{ ^ . (139) 
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Because E'f(r) has its source in fF^) and this is excited by E^f"), the interference is 
coherent. 

In experiments, it is possible in principle to measure both the total and the coherent 
intensity. One can think for instance of a medium made of dipole antennas in which one 
of them is excited by some external fixed field of frequency u = ck . Part of the power 
supplied is absorbed directly by the source antenna located at r, which is given by the term 
of Eq.(41). The rest is radiated. Again, part of the radiation will be extinguished, but 
W Coh -(r) as given in Eq.(131) will propagate coherently through the medium and can be 
collected, eventually, by some distant receiver. If the medium can be treated as an effective 
medium, the coherent intensity collected by a receiver located at r' will be given by 

W Coh -(\r' -r\) = W Coh -(f) exp [-2k \r' - r\R], (140) 

where R = ^{\/Xe//} is the extinction coefficient which contains both dispersion and ab- 
sorbtion by the host scatterers. 

It is also instructive to compute W Coh - for the stimulated emission of one of the atoms of a 
monoatomic dielectric close to the resonance. Let us assume the polarizability of the atoms 
adjusts to the Lorentzian function of Eq.(44). Using Eq.(40), we can write the coherent 
power emitted by the stimulated atom sited at r as 

W at°l = -^l«W| 2 IK| 2 Tr-{ / dV»{*^}3{G(*V;a;)}}| , (141) 

DC p <- J Ct J \k=u)/c 

Close to the resonance, a(k res ) ~ i ck res a /T is nearly pure imaginary and the total power 
supplied reads 

= ^y^\E^\ 2 Q{&(k res )} w ^p|£ — |2|d|. (142) 

We can write < 3{G(r',r;u)} in terms of using the fluctuation-dissipation relation of 
Eq.(132) and, making use of the constitutive relation for coherent fields, write the polariza- 
tion vector as 

P"(f) = e J dV x(r,r) ■ E%{?). (143) 
Finally, by inserting the above results into Eq.(141) we obtain the more familiar expression 

WZt = -^P£%{(-P** T "(f) ■ K-tm, (144) 

— * 

where i-P* w ^{f) is the average molecular dipole moment. Therefore, we see that in the VC 
scenario there is a well defined ratio between the coherent emission of a stimulated dipole 
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close to resonance and the total power supplied. That ratio being given as a function of 
coherent fields. 



2. The coherent vacuum and the additional polarization of the sourceless vacuum 

In the following we interpret the coherent emission in QFT terms. From Eqs.(132,136), we 
conclude that the states accessible to direct coherent emission are the normal modes which 
satisfy Eqs. (125,126). Therefore, in an effective manner, the vacuum though which these 
modes propagate is the sourceless vacuum, |fi) sJ \ This effective equivalence is analogous to 
that found in Section II between \£l) s {n and |0) for a unique dipole in free space. \Vt) s - L is 
made of the superposition of EM vacua, {\(j>m) s ' 1 '}, in which the normal modes of Maxwell's 
equations with permitivities {€ m (f)} propagate. Following the formalism of Section II C, 
the m th pair of vacuum state \4> m ) s ' L and dielectric permitivity e m (r) is associated in a 
one-to-one correspondence to the m th configuration of host scatterers such that \il) s - L = 

E m v / M^:i0 m ) s -s 



J2 s - L {<t>rn\M^{^{r)\<t> n ) s - 1 - = s - L (Q\^(f)^\r)\n) 



s.l. 



= -h^^{G(r,7-u)}. (145) 

6q7TC 



On the other hand, it is plain from Eqs. (127,129) that the modes which contribute to 

W Coh. 

are the same as those which propagate in |fi) sX . However, the amplitude of the 
spectrum of fluctuations which enter W Coh - differs w.r.t. that of the fluctuations in |f2) sJ - 
by a multiplicative factor ^{x/poi}. In other words, the EM vacuum in which coherent 
emission propagates, \Q) Coh , is additionally polarized w.r.t. The renormalization 

function is given by y = ^|-^~"-|, which takes account of the polarization due to the 
closest scatters surrounding the emitter. We will see later on that y is the so-called local 
field factor in the theory of the effective medium. In the framework of QFT, y can be 
interpreted as the field-strength renormalization factors. In a scalar QFT it is the analytical 
structure of G(k) together with Z in the complex plane that determine the amplitude of 
the vacuum fluctuations and the mass spectrum of the particles propagating in space-time 
[29]. Correspondingly, it is the analytical structure of G± j \\(k) together with 3t{ X± ^^ } that 
determine the amplitude of the EM vacuum fluctuations and the /c^n^'-spectrum of the 
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coherent photons propagating in a random medium. Following up the QFT interpretation, 
^± || W yields the probability density for creating a photon of frequency u> from vacuum at 
the emitter site, Z u) (k) = \ Coh (Q\E ± ^(r)\^)\ 2 . Out of those photons, the coherent ones are 
selected and propagated through the bulk by G°[. Also, the electric field gets renormalized 

~* —1/2 ~* 

and the field so renormalized is nothing but the Dyson field, I_l,|| = Z± ^ E^_ y. This leads 
to the equivalence relation, 

sJ -(fi|l"(f)l wt (r)|fi) si - = Coh {n\^(r)^(f)\n) Coh 

= -h-^Z{G(r,f<-u,)}. (146) 

— * — * 

Thus, in passing from E% to E u , any correlation function between vacuum states picks up 
a factor [Z"] 1 / 2 per field operator. 

All this suggests that an effective electromagnetic theory in random media can be formu- 
lated in terms of the renormalized quantities -see eg. Ch.7-10 of [29]. The host scatterers 
play the role of bare vertices of photon interactions as depicted in the diagrams of Fig. 15. 
The correlations, depicted by dashed lines in our diagrams, play an analogous role to that 
of the gauge fields which mediate interactions between charged fields in QFTs. If corre- 
lations are assumed not to have dynamics in accordance to the weak-coupling-Markovian 
approximation, they must be intended as a back-ground field. 



3. Radiation in an effective medium 

In this section we deal with the effective medium theory, which is implicitly used in the 
majority of works on dipole emission. In this context, and only in this context, an analog 
contribution to the radiative emission in free space can be defined. 

The effective medium theory is defined such that, for some range of frequencies and for 
fc( < 1, an effective homogeneous complex dielectric constant can be defined as e e // = 
Lim. {e±,||(fc)} as k£ — > 0, so that the Dyson propagator G±(k) can be approximated by the 
effective propagator G e /*(k) for some range of frequencies, 

Gl ff (k) - - * , Gf{k) = * (147) 
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The field with propagator G e // will be referred to as Dyson-effective field or macroscopic 
field, . Correspondingly, the field with propagators 



pa e eff k^-k 2 " pa e eff k% 

will be referred to as local-effective field EfJ/. Further on, in order for G e /*(k) to be the 
Green function of a propagating field of wave number k it is required k/k > 1 at the same 
time that ^ < 1, which implies k ^ C 1 as well. 

To what computation of power emission is concerned, the use of macroscopic fields implies 
the neglect in the integrals of Eqs. (127-130) of the poles of order > l/£. The neglected 
modes correspond to the near field dispersion and absorbtion carried out by the scatterers 
surrounding the source -see Fig. 16 and also [15]. In turn, this implies the neglect of all 
longitudinal emission and the restriction of the integration domain in Eq.(123) for transverse 
emission, 

WJ> S 2^- * W.J*^ *{*%-^} (149, 

= 2W^{^±^}, W^~W C . (150) 

The coherent component of the above formula, 2W^*\ Coh - can be written in an analogous 
manner to the radiative emission in free space, 

2W eJffoh. = 2HMftj^ js R { v /^}. ( 151 ) 

2W e ± ff \ Coh - can be given in closed form for a Maxwell-Garnett (MG) fluid. In such a fluid 
the only relevant correlation between host particles is the negative correlation due to the 
exclusion volume ~ (47r/3)£ 3 which prevents scatterers from overlap. Provided k Q ^ <C 1, 
e e ff obeys the MG formula [58], pa = ^ €e j f f f +2 y 3 - It was firstly proved in [15] and then in 
[59] following a diagrammatic approach that the MG formula is correct at all orders in pa. 
Note however that the usage of the renormalized polarizability a incorporates in our case 
the effects of recurrent scattering. Using the MG formula, Eq.(150) reads 

2W M G] Coh. = 2^| £e ^ +2 |^{ v ^7}, (152) 

It is in order to make a couple of comments regarding the use of effective fields in the 
computation of LDOS, W and V. First, in Eq.(149) the symbol ~ stands for the approxi- 
mation SoKhKA-H- 1 (£f ~ fk>o (£f under the assum P tion ( fc oO _1 > A > L That implies 
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that the integral of Eq.(149), which is the inverse Fourier transform of G± (k) evaluated at 
r = 0, does not actually carry information about the fluctuations of the microscopical local 
self-polarization field evaluated around the emitter location, r < £. Instead, the fluctuations 
there are those of the field Ef„/ , which is the local field spatially-averaged in volumes of 
the order of £ 3 over which also the spatial average of e(r) is performed in order to get e e //. 
Therefore, because the contribution to W\\ comes from high frequency modes in the range 
k ~ V£> ft i s n °t possible to take the effective dielectric constant limit k£ <C 1 for its 
evaluation and it is in general unavoidable to incorporate the spatial dispersion in x\\(k)- 
Differently to normal propagating modes, the electric field of the source dipole is not ap- 
proximately uniform within a cavity volume of the order of ~ £ 3 and the emitted photons 
do not see an effective (uniform) dielectric constant there. In the real large-cavity scenario 
however, it might still happen that the field at distances 1/R is approximately the same as 
that at 1/(R + £), £ being the typical dispersion length scale of x m the host medium. It is 
only in that case and for fc £ <C 1 that it is possible to consider x uniform outside the real 
cavity. 

Second, let us assume for a moment that the effective field approximation is good enough. 
It is known long ago that the effective field which appears in Fermi's Golden rule in the form 
of the equal-point correlation functional, (E(r)E(r)) , is not i?^ (r) but E^/(r). Further 
on, assuming the medium is passive, emission restricts to 2W^^\ Coh \ In application of the 
fluctuation-dissipation theorem, 

2W e ± fffoh. = W Q j 0^^{^}^{ G ± f ( k )} 

= ^f^ { ^ }( ^ (r) ^ /(r)) - (153) 

Notice that it contains just a factor ^{^-f}, which equals a Lorentz-Lorentz local field 
factor for the case that MG formula applies, C LL = ^{^f^}. It is a common mistake to 
think that Cll should appear as squared in front of {E e p \r)E e ^ (r)) in Fermi's Golden rule. 
The reason being that while the propagator of E^/(r) is that of E e J^(r) multiplied by £ LL , 
Q e ff ~ CihGeff, the fluctuation-dissipation theorem applied over Q e ff yields a contribution 

— * 

quadratic in E D but linear in Cll- This contradicts the paradigm which states that the 
effect of a dielectric host is to multiply each occurrence of the dipole moment of a a two-level 
atom by a local-field factor [22] -see Section IX for further discussion. 
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4- On the radiative and non-radiative energy transfer 

Taking advantage of the above study, we can interpret the radiative and non-radiative 
nature of the processes involved in the dipole emission of the single scattering model. Let 
us consider the diagrams in Fig. 16. They represent either the stimulated emission or the 
spontaneous emission of an emitter in the single-scattering approximation -Eq.(103). In the 
case of stimulated emission, /i 2 stands for the polarizability assigned to the emitter. The 
propagators on the right hand side of the emitter are part of the bulk propagator G. Actual 
emission is denoted by wavy lines while induction is mediated by propagators depicted by 
solid lines. Curved propagators on the left hand side of the emitter indicate being con- 
strained by spatial correlations within x- Emission with origin at the emitter is direct. It is 
induced emission otherwise. 

Since the work of Andrews [28, 60], it is well appreciated that the radiative and non- 
radiative energy transfer between a donor and an acceptor correspond to a unique quantum 
mechanical process. When no induction exists, the emission is direct and there is no trans- 
fer of energy. It corresponds to the diagrams (a), for in- free-space emission; and (6), which 
stands for single scattering in the bulk -Eq.(113). 

Induced emission processes, either coherent or incoherent, are characterized by two fea- 
tures. First, there exists some correlation between the donor and the acceptor. Second, it 
is one of the transverse virtual photons of frequency u res which mediated the interaction 
donor-acceptor, that is made actual when satisfying the dispersion relation of Eq.(125). The 
pair donor-acceptor behaves as an effective dipole. That in diagram (c) is coherent and so 
is its (transverse) radiation, as correlations apply to the electrostatic induction -Eq.(115). 
On the contrary, that in diagram (d) is incoherent, as correlations apply over the transverse 
induction whose photon is the one becoming actual -Eq.(118). As explained above, the 
processes in (c), (d) carry one local field factor each in the single scattering approximation. 
They are the processes which can be considered to amount to radiative transfer. 

The non-radiative transfer processes are those in diagrams (e, /, g). They are character- 
ized by two features. First, the actual photons have their origin either in radiative correc- 
tions over the acceptor bare polarizability or on the imaginary part of the bare polarizability. 
Hence, the corresponding terms are all proportional to $5(0;}. Only radiative corrections are 
depicted in Fig. 16. Second, the transference of energy takes place through either purely 
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-diagrams in (e)- or partially -diagrams (/, g)- longitudinal loops which remain closed 
during the re-emission process. In general the spectrum of energy levels of the acceptor 
presents a number of channels into which the quantum of energy absorbed, hjj res , can be 
re-emitted. Therefore, re-emission takes place in a variety of frequencies and the associated 
decay rate is said non-radiative. Note however that this implies a slight inconsistency in 
the case that donor and acceptor be identical two-level atoms. That is so because in such 
a case the acceptor can only re-emit in the same frequency as that of the direct radiation. 
Nevertheless, because two-level dipoles are unrealistic idealizations, we will consider those 
processes as non-radiative without loss of generality. The diagrams in (e) stand for the 
terms in Eqs.(116,117). They would diverge if considered separately. They represent the 
phenomenon commonly referred to as Foster Resonant Energy Transfer (FRET) [61]. The 
diagrams (/) and (g) correspond to the absorptive longitudinal and transverse terms of 
Eq.(118) and Eq.(115) respectively. 

VII. ANALYTICAL CALCULATIONS IN MAXWELL-GARNETT DI- 
ELECTRICS 

An MG dielectric is one made of point dipoles with well defined single particle polarizabili- 
ties. The correlation length £ among dipoles satisfies fc £ <C 1 for the frequencies of interest. 
Under these conditions, the only relevant correlation function between its constituents is 
given by an exclusion volume correlation function. In such a medium, an effective dielectric 
constant can be defined and adjusts to the MG formula, 



Hereafter we will denote e e // with 6mg instead. In this scenario it is possible to obtain ana- 
lytical formulae for the virtual cavity 7-factors. We apply those formulae to two situations of 
experimental interest. In the first one, we compute the decay rate of an excited atom which 
can be either an interstitial optical center or one of the molecules of the dielectric itself. In 
the second situation, we compute the effective dielectric constant for an MG dielectric close 
to the resonance frequency. 
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FIG. 16: (a) Diagrammatic representation of the stimulated emission/spontaneous emission of an 
emitter in the single-scattering approximation -Eq.(103). Explanation is given within the text. 

A. Spontaneous decay rate of a weakly-polarizable excited atom 



Let us consider the explicit computation of T for an excited atom. This can be either a 
foreign ion placed within one of the molecules of the dielectric or a molecule of the dielectric 
itself. In the former case, the manner the emission rate of an optic center gets modified when 
placed within the ligand is a difficult issue as, in general, one should consider the internal 
atomic structure of the host molecule [62-64]. To simplify matters, we assume here that the 
bare dipole transition amplitude of the ion and its resonance frequency get renormalized to 
known values result of the interactions with other atoms within the ligand. This 

scenario is that of Section III C. The electrostatic polarizability of the doped molecule, a , 
is carried out by some other polarizable atoms, being a equivalent to the bare polarizability 
of the rest of the host scatterers. The total dipole moment associated to the particle which 
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2 — 

hosts the ion is the combination of a fixed and an induced dipole, phost = p + ^■aQ UJ {f,f)p, 
a being the renormalized polarizbility of the host particle. 

The case of an excited host particle can be treated in the same fashion provided that 
the additional condition £ <C p -1 / 3 is satisfied and the effective polarizability of the excited 
molecule is negligible -see Section VA. That corresponds to the setup of Section III C with 
a = and assuming previous knowledge of p and uj res . 

In this scenario, the only correlation function which matters is the two-point function 
which amounts to the exclusion volume of each atom, h(r) = — 0(r — £). For higher corre- 
lation orders and following [15] we will make use of the overlap approximation. Ultimately, 
because Eq.(154) establishes a relationship between xmg and a, and the 7-factors can be 
expressed as power series of pa, a relationship can be given between T, LDOS and the ef- 
fective parameters which determine the transport of coherent light. Those are the index of 
refraction, n, and the extinction mean free path, l ext , which relate to Xmg through 



n = ^Wl + X M G }, l; x \ = 2^{y/l + X MG}- (155) 
First, we compute Q±. In the overlap approximation, 

X ^ = \ WWTi -v (156) 

l-xT(k)/(pa) 

where at resonance, 

X { l\k) 



(2>, M . MXes I d 3 re ^( r ) Tr {G(0)( r )[I_ k(g) k}} 



2 



h(\k'-k\) Gf\k') 



2 J (2vr) 3 

+ (k r ) cos 2 9 + G§ ] (k f ) sin 2 0~\ . (157) 

For an MG dielectric, h(Q) = -^ 3 ji{Q)/Q = 47r£ 3 Qcos g 3 ~ sinQ with Q = \k' -k\£ and 3x the 
spherical Bessel function of first order. Applying next Eq.(58) at the resonance frequency, 
we find 

Q ± (k) = Gf\k) \l + pa[kl es Gf{k) - xf{k)/{paf]} \ (158) 



From the above expression it is not straightforward to differentiate between the coherent 
and the incoherent components of 2r_i_ as given in Eqs. (127,128) for the power emission. 
However, it is possible to go around this problem by first computing Eq.(158) order by order 
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in ( = k res C, and then find out the origin of those terms in view of Eqs. (127,128). First, we 
write 2 7 J°*- as 



7 POO f~\1 00 /~1 

2tIo , = ^| ^ dg _^£ ( _^ r[ _i_ - (pa)- 2 xi 2) (CO)]"}, (159) 

n=0 



vr 2 CU "C 2 -Q 2 ±f r y L C 2 -Q 



where Q = k£. Careful examination of the integrand of Eq.(159) yields the following de- 
composition which includes all the terms up to O(0) in (, 

27f( . K ^^.^ fy Q ^_ { ^ _ 1/3] . (160) 

' n=0 u ^ 

+ E(-^) n / d ^b'i(Q)/Q] n } (i6i) 



n=l 

= 2 7 1°'- A + 2 7 l '- B . (162) 

The integrand of Eq.(160) contains the poles of G± which satisfy the dispersion relations of 
coherent modes. Eq.(161) carries those very high frequency modes, k ^> k res , which contain 
the near field contribution. By taking the limit ( — > 0, it is immediate to identify in Eq.(160) 
the factor —1/3 (pa) 2 with the zero mode of — x±\k). In the overlap approximation, that 
leads to the MG formula for the effective susceptibility. Hence, it is x^i^ = 0) = 1/3 (pa) 2 
that enters in the place of x±\k) in Eq.(158) for Q e /^ . We take advantage of this to arrange 
the leading order terms of Eq.(160) (which are those of O(0) in (). We obtain, 

2T A ± = 2f 3{-2 7 r^ « 2f ^{^±^ v ^), (163) 



v res 



k 3 



^-k 3 a + h 2 a [2-f ± + ^ (164) 



where T Q contains the 7-dependent factor which renormalizes \p\ 2 due to the polarization 
of the host particle. It is also immediate to verify that the coherent contribution of 2Tj_ is 
precisely 2T^ G \ Coh - as given in Eq.(152) (modulo prefactors), 

2T A,coh. = 2V M G] Coh. = 2f Sft{ eMG 3 +2 }sft{ v /^}, (1 6 5) 

while 



AJncoh. _ -p o / e MG + 2 



-}%{Vq^}. (166) 



°*Y 3 

Therefore, at leading order in ( = k res C,, we conclude that the coherent emission -and only! 
the coherent emission- can be computed in the framework of the effective dielectric constant 
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approximation. 

On the contrary, incoherent emission does not only contains the term 2T^_' Incoh ", but also 
the terms in 2T I [ which are proportional to Q{2'y-[° t - B } given in Eq.(161). 2 r y 1 [ ot - B takes 
account of those non-radiative high frequency modes which scale as Eq.(161) contains 
terms of the form 



"OO 



7 $ = / d Q 1^ with P < 2? + 1, 2n < p < 3n, (167) 



'o 



J&> = jf dQ "-^1 with p < 2q, 2n < p < 3n. (168) 

With increasing value of n ^> 1 it is easy to verify that the above integrals converge quickly 
to zero. Hence, 

/« = n - { - 1)9p2q (169) 
2 (2g)! 1 j 

J« = iW**- 1 . (170) 
p ' 9 2 (2g-l)! { } 

Therefore, for any finite value of pa, the corresponding series in equation Eq.(161) is rapidly 

convergent. We give below the first five terms of that series, 



—k 1 
2^ B ~ ' '" " ' 



1 _ 2 . . 2 47 . _, 3 334 . , 4 6891623 . _, 5 
2^-15^ + 1^^ -31185 ^ ) + 214^927168 M + '-- ' ^ 



27T ( 

Let us consider next the computation of 7y. The Green function for longitudinal modes 
reads 

^||(A:) = T^^[l + Xi|(A:)]- 1 , (172) 

^res F u 

where in the overlap approximation, 

xnW = - — m^jr- v ( 173 ) 

1 - X {k)/{pa) 



with 



x j 2 )(A;) = -(ap) 2 k 2 J &\e it %r)Ti{G {a \r)[k®k}} 

J 3 7. / 



X 



G^^Ocos 2 ^ + Gf(fc')sin 2 . (174) 



The above equation can be computed in closed form, 

X f{Q) = {paf[l + ^/(C)], /(C) = 2. e*(i + 0- (175) 
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Likewise, it is possible to give a closed expression for G\\(Q) using the MG relationship for 
€mg m function of pa, 

ai(Q) = ^ 1 + e 7 !£ i^^' 1 - (176) 

«Ves e MG + z W 

The use of this expression in the computation of T\\ implies the computation of an infinite 
series of roots of 1 + 6 emg+2 ^P -'■ ns ^ eac ^ °f that, we rather expand Eq.(172) in powers 
of pa as this allows to keep control over the convergence of the series, 



W) = T^J2(p & n-i) n no[h(Q)/Q\ n - (m 

Finally, we arrive at 



k 2 

res n=0 



7|| 



k 



1 _ _ r°° 

> n=l ^ 



Note the similarity between the above equation and that for 2'yJ_ otB in Eq.(161). This is at 
the root of the coincidence of the terms of order pa/( in the single-scattering model studied 
previously. The convergence the series in Eq.(178) can be verified using similar arguments 
to those employed for Eq.(161). 7y contains however additional terms of order and 
order zero in ( as a result of exponenzing the function /(£) in Eq.(175). We give below the 
first five terms of the orders £>(C" 3 ), ^(C -1 ) and £>(C°), 

* + 1^ + + + + -1 < 179 > 

k res lA _ 2. 5, 544. _, 4 204745, ... . 

— t -P« + -(pa) 2 + —(pa) 3 + (pa) 4 + (pa) 5 + ... (180) 

2tt (2 3^ ; 16^ ; 2835^ ; 1741824^ ; 

.A; res 1 _ 4 2 5 3 1088 4 204745 5 



x- 



pa + -(pa) 2 H (pa) 3 H (pa) 4 H (pa) 5 + ...1. (181) 

2tt l 3^ 9 V ^ ; 24^ ; 8505^ ; 2612736 yH J 

For an MG dielectric, it is possible to distinguish between dispersion and absorbtion in the 
same manner we did for the simple single-scattering model in Section VI B 4. In particular, 
the terms in Eq.(179) are the analog of a FRET process for the case that the transfer of 
energy takes place between the emitter and an MG dielectric. Below, we write the complete 
formulae for the different contributions to the decay rate (up to renormalization prefactors) 
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in function of the effective dielectric constant, 

r£& = 2 f ^{^|±^{ v /^}, (182) 



r lDis P er. = 2 f Jk{ ]_ XMG + I X ^ G _ o.051x^ G + 0.055 xl, G - 0.015 x^ G + -..} (183) 



-pAbsorb. 
1 MG 



= -2f S{^^}s{v^} (184) 
of 1 1 

+ ^{ 2^mg + - 2 X 2 mg ~ 0-076 x^g + 0.083 X \ig ~ 0.022 X ^ G + ...} (185) 

2f 11 

+ ^{Xmg + -Xmg - gXMG + 0-073 Xmg " 0-028 + •••} (186) 

Note however that when multiple-scattering processes are involved, it might well be that 
part of the decay attributed to absorbtion is actually due to dispersion by clusters of cor- 
related host dipoles. A microscopical analysis of the origin of ^s{xmg} would be needed to 
discriminate between both extinction processes. 

At this point, it is interesting to compare the order O(( ) term of T^Jf to that usually 
employed in the literature which includes the square of the Lorentz-Lorenz local field factor, 
T LL = 2f (^§±2) 2 v /i^G". In our computation, = T c ^ + Y^ sper - . These terms are 
the relevant ones for emission in absence of absorbtion, ^s{eMG} ^ ^-{^mg}- Expanding in 
powers of Xmg — € mg — 1, we obtain up to order five, 

7 3 



C G = 2f o K{l + - XMG + - X ^ G - 0.030xLG + 0.037x 4 MG- 0.0007 X ^G + -}(187) 
T LL = 2f o ^{l + ^XMG + §XMG + 0.035x 3 MG-0.01xMG + 0.008xl f G + -}.(188) 



Therefore, in absence of extinction, the agreement is almost perfect up to order Xmg- 

As 

we will argue in the next section, this agreement is just accidental. Therefore, in order 
to differentiate between both formulae in non-absorptive media, it is necessary that higher 
order terms be relevant. That is, the index of refraction has to satisfy n > y/2. Never- 
theless, we must point out that it has been assumed implicitly in Eqs. (187,188), by using 
a unique value for r o , that the renormalization factor which accompanies fi is the same in 
both equations. For the decay of a weakly-polarizable host molecule this is not a problem. 
However, for the decay of a foreign particle within a host molecule this might not be so. 
Even for a non-dissipative medium the terms of orders 0((~ 3 ) and 0(C,~ l ) are relevant (see 
Eq.(49)) and they do not appear in the Lorentz-Lorenz approach. Therefore, any try to 
differentiate between both approaches must necessarily assume some ad hoc common value 
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for the renormalization factor. We give below the 7-factors which enter the renormalization 
factor in Eq.(164) bearing in mind that such a renormalization factor must be set to 1 for 
the case of a weakly-polarizable excited host molecule, 

27_l + 7|| = 7c° + 7c- 1 + 7c- 3 (189) 
k T 3 

= -^qXmg + ^Xmg - 0-030 Xmg + 0-037 X mg ~ 0-0007 XmG + -1(190) 

- ^[[Xmg + \xmg + 0.105 xl G ~ 0-027 Xmg + 0-006 XmG + •••] (191) 

- ^^Ixmg + IxIg - \x 3 mg + Xmg Xmg + ■■■}■ (192) 

The above formulae are thought to fit the experimental results analyzed quantitatively in 
[65]. There, Ce +3 ions of low polarizability replace low-polarizability cations of several host 
media. The life-time of the transition 5d — > 4/ is measured in several hosts and the data 
are fitted in [64] to the usual VC and RC formulae. The authors find 'relative' agreement 
with the VC model. However, in their fits they leave the value of pi to be fitted as well, 
introducing this way an additional degree of freedom which improves artificially the fits. 
Therefore, their results are far from conclusive. 



B. The dielectric constant of a Maxwell- Garnett dielectric 

The problem in this case is one of self-consistency. On the one hand, we must compute 
the renormalized single-particle polarizability of the molecules, a, due to their embedding in 
the dielectric. On the other hand, the renormalization of a gives rise to a renormalization of 
the dielectric constant itself with respect to the MG formula involving in-free-space polariz- 
abilities. Following the renormalization scheme developed in Section III B, beside the value 
of «o an d the in- free-space resonance wave number, fc , knowledge of the exclusion volume 
(i.e. the correlation length £) and of non-radiative effects (i.e. the collisional shift A£^ oH 
and collisional line broadening T coll ) are needed. While the latter effect are thought to be 
relevant at high temperature, radiative effects are expected to dominate at low temperature 
and for frequencies close to the resonance. 

A possible setup which adjusts to this scenario corresponds to that in the experiment 
carried out in [66]. There, selective-reflection techniques are employed to measure the fre- 
quency shift in a high temperature potassium gas. At the frequencies of interest, the bare 
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electrostatic polarizability reads a — f Anr e kQ 2 where r e = 4 ^ c2 = 2.82 10~ 6 nm is the 
electron radius, the resonance wavelength is A = 2nk 1 = 770. lnm and the strength factor 
is / = 0.339 for the transition 4 2 Si/ 2 ^ 4 2 Pi/ 2 [67]. The range of atomic densities in 
the experiment is 10 23 m -3 < p < 2 • 10 20 m~ 3 . We can keep the renormalization scheme of 
Section III B by just replacing a with a renormalized static polarizability which includes 
both the regularization of 29?{ 7 |°" ) } and collision factors, 



<°tat.(k) = a k 2 [k 2 + Ak 2 colL - i T coll ck - k] . (193) 
This leads to the Lorentzian renormalized single-particle polarizability, 

a(k) = a k 2 [A; 2 + Ak 2 olL - i T colL c~k - k 2 + ±a l4k?(i$t{2>y ( ®} + 2 7± + 7,^] \ (194) 
Next, in application of MG formula, the effective dielectric constant reads 

e MG (k) = 1 + - Pa } k } r ~ s = 1 + XMcCk) = 1 (195) 



1 — irpa(k) 



+ pa Q kl ]kl + Ak 2 coll - i T coU k/c~ k 2 + ^a klk 2 (i^s{2-f < f ) } + 27 ± + 7^ - ^a k 2 p 

The 7 -factors are those of Eqs. (189-192). Because they depend on Xmg, the problem be- 
comes one of self-consistency. Self-consistency can be seen also as a consequence of including 
recurrent scattering in the derivation of the MG formula. In physical grounds, it reflects the 
complementarity of the double role played by the actual dipoles. That is, on the one hand 
they polarize the sourceless EM vacuum characterized by cmg- On the other hand, they 
renormalize their own polarizability, a. 

Once the self-consistency problem has been solved in terms of Xmg, the 7 -factors can 
be computed and the renormalized values of T, oj res and ao follow the equations of Section 
III B. Note that the resonance frequency of l MG differs w.r.t. to that of a in the so-called 
Lorentz-shift, Ak\ = —^a klp. 

The authors of [66] however did not include all the possible radiative effects in their anal- 
ysis. They restricted themselves to the Lorentz-shift, Ak\. Because Ak\ does not depend 
on Xmg, no self-consistency was demanded. In the experimental setup of [66] no further 
radiative effects are relevant despite the fact that the medium is highly opaque. The reason 
being that collisional effects dominate. A detailed analysis of the problem -see [68]- reveals 
that this is the case. The main reason is that those ( dependent terms in Eqs.(191,192), 
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which would be expected to become large for \xmg\ ~ 1, *C 1, turn out to be irrelevant 
in comparison to the collisional contributions. For a cold gas of potassium, a estimate for 
£ would be the van-der-Waals radius of K atoms, which is roughly 0.5nm. However, at 
high temperature £ is rather determined by the collision cross section between atoms, which 
turns out to be much greater. It would still be possible to increase the relative weight of the 
radiative terms with respect to the collisional ones by lowering the atomic density. However, 
for a precise calculation we would still need a good estimation for £. 

To the opposite extreme belongs the scenario in which the medium is an atomic glass. In 
such a case, there is the advantage that £ can be accurately determined. In fact, £ ~ p -1 / 3 . 
However, as pointed out in Section II C, atomic orbitals overlap in such a way that the 
electronic band configurations can be very different to those of the individual atoms iso- 
lated. On top of that, there might be contributions of free electrons, which introduce an 
additional source of uncertainty in computing that component of the susceptibility due to 
plasma oscillations [69]. 

In view of the difficulties that high temperature atomic gases and solids present in evalu- 
ating the radiative contributions to the renormalization of the dielectric constant, we propose 
alternative setups. They must be such that undesired dynamical and collisional effects be 
reduced to a minimum. In addition, the electronic structure of the dip ole- constituents must 
be hardly altered by their cohesive interactions. A cold atomic gas offers a scenario which 
adjusts to these requirements [70]. The other possibility is the preparation of colloidal liq- 
uids [71]. There, particles of well defined polarizability are in suspension on a background 
solvent. The particles are ionized in such a way that the spatial correlation among them 
is determined by the Coulombian long-ranged mutual repulsion and crystallization is sup- 
pressed. The degree of order and the value of £ can be tuned to satisfy the conditions of an 
MG dielectric. 



VIII. SPONTANEOUS DECAY OF A LOW-POLARIZABILITY SUBSTANTIAL 
IMPURITY 

In this section we compute the decay rate, T, of a fluorescent weakly-polarizable emitter 
with transition dipole moment \x and resonance frequency u res = cko. The emitter locates 
at the center of a large cavity of radius R much larger than the correlation length £ among 
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(a) 



FIG. 17: (a) Cermet topology scenario. The emitter with polarizability ao is surrounded by a 
disconnected net of scatterers with correlation length £. (6) Simply-connected-non-contractible 
topology in which the emitter is placed at the center of a cavity of radius R surrounded by a 
continuous medium of dielectric constant e. 

host scatterers. The polarizability of the emitter is low in comparison to that of the host 
scatterers. The long- wavelength limit applies on the host medium, /c £ <C 1, such that it 
behaves locally as an effective medium of uniform dielectric constant e = 1 + x- The setup 
corresponds to that of Section VB. The topology associated to the embedding of the emitter 
is that of a simply-connected-non-contractible manifold -see Fig. 17(6). Its experimental 
counter part corresponds to that of Eu 3+ ions in random media [72] . Our master formula is 



where the 7-factors are those of the RC scenario given by Eqs. (81-96). Following those 
equations, we will compute the 7-factors as power series of x- Our computation is exact 
up to order two in x- 111 order to compare with the conventional empty-cavity model 
formulation, we will restrict ourselves to the case k$R <C 1, although exact formulae at all 
the orders in k^R can be obtained beyond this approximation. 



r = 



3c 2 h 



H 2 3{2 7 r + 7 f c }, 



(196) 
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The 7-factors are given by 

+ -^- 3 [C ± + G^rG ± xl(k), (197) 



+ k 4 J ^[C\\+G\^G^(k), (198) 

where the cavity factors are those of Eqs. (90,91). Considering that the polarizability of the 
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FIG. 18: (a) Feynman's rules . (6) Diagrammatic representation of the transverse and longitudinal 
7-factors, and 7jj ?C 'j as given by Eqs. (197,198). 

emitter is negligible and that the interaction of the emitter with the surrounding particles 
does not modify its energy levels, neither changes on /i nor in the resonance frequency hold. 
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That way, we can perform the substitution of k with ko in the above formulae. Up to zero 
order in k R we obtain, 

-e + 2\2 n 4(e-l) 2 " 



2 1± C + 7jf ° - 



k 



2tt 

ko_ r 

2tt 



9 



1 



+ 
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e + 2, 



.(M?) 3 k$RJ LV ~ y 3e 
Normalizing by the in-free-space decay rate, T = 3 ^ e e 3 3;i |/x| 2 , V reads 

'e + 2\2 
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(201) 



.(Jfeo-R) 3 ' A; 0j R/ " L 3e 
The terms in Eq.(201) are associated to absorbtion in the host medium [73]. We recognize 
in the first term of the r.h.s of Eq.(200) the usual bulk term corrected by Lorentz-Lorenz 
(LL) local field factors, T LL = T ^{(^4^j y/e}. The second term there includes corrections 
of order > \ 2 ■ It is also remarkable that if the empty-cavity Onsager-Bottcher (OB) local 
field factors are used instead, T^f = ro3^{^2efi) V^} does agree with the two radiative 
terms of Eq.(200) up to order x 2 , 

7 1 



^emp 
OB 
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-X + ■■■)■ 



(202) 
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Our computation is in fact exact up to 0(2). However, the nature of the emission associated 
to the terms of Eq.(200) is different to the one attributed in the conventional empty-cavity 
model. In the latter, is considered as fully transverse. We show next that transverse- 
coherent emission is just part of it, being the rest incoherent longitudinal emission associated 
to dispersion. We can read the contribution of coherent emission from the diagrams in Fig. 18, 



^Coh. 
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(203) 



where the first term on the r.h.s of Eq.(203) corresponds to direct-coherent emission - 
transverse components of diagrams (bi), (b 2 ), (64) in Fig.18. In the above expressions, 
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C±(y/ek ) stands for C±(k)\ k=y ^ ko . Y Coh - can be written in a more compact form as 

„ -2tt 



^Coh. 



— [2»{[1 - kl X C ± (V~eko)} 2 } J %{G^)}-^ (204) 
+ 2K{k 2 oX C ± (V~eko)} J Q{G^\k)}-^\, (205) 



which immediately leads to 



j^Coh. -p 

1 RC — 1 



nv~e}®{ (4^) 2 } - ^ e - x >] > ( 206 ) 



at leading order in k R. For the dispersive and absorptive terms we get, 



r 



,Disper. 
RC 



.*{^§^}. (207, 



rgT* = -r„[3{v^}3{(^-) } (208) 

+ ((W^) 31 ^-*'- 11 '] (209) 

Once more, the same reserves after Eq.(185) about the distinction between absorbtion and 
dispersion apply. 



IX. COMPARISON WITH PREVIOUS APPROACHES 

We have proved that our computation of LD 0S^™ SSJori of Eq.(32) with the 7-factors 
given by Eqs. (60,61) is exact for a microscopically linear and statistically homogeneous 
dielectric in which the electrical susceptibility tensor can be expanded as in Eq.(55). Here 
we review on previous works, compare them with ours and emphasize where the errors in 
those approaches are. Because the literature about complex media concentrates in setups 
where an effective medium can be defined, we will restrict ourselves to such a particular 
scenario. Our analytical study of an MG dielectric of Section VII belongs to that kind. 

A. The erroneous inclusion of local field factors 

The majority of approaches on non-dissipative media agree in giving an effective density 
of states of the form 

LD0S erms Sl on = |£(^ oc .(f) • E$J?)) » \ k<< ^C) U LDOS° ut (210) 
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where Cf is a local field factor. The approximation symbol denotes the restriction to long- 
wavelength modes with k <C which is implicit in the use of the effective medium 
theory. (E^ oc (r)) is interpreted as the electric field felt by the emitter located at r and 
(^Loc if)^i}oc (0) * s ^ ne s P ec trum of its vacuum fluctuations. As pointed out in Section 
VIB3, because modes with k > are not included in Eq.(210), such a field is actually 
not microscopic but assumed uniform within the cavity of radius £. Hereafter we use the 
nomenclature of Section VI B 3 where E\H \r) stands for E^ oc (f) in the effective medium 
theory. Likewise, E^^(r) denotes the macroscopic field in the bulk. The usual procedure 
in the literature is to compute E^/(f) classically by imposing boundary conditions at the 
cavity surface. Those boundary conditions depend on whether the medium inside the cavity 
is equivalent or not to the medium in the bulk. In both cases, effective media are assumed to 
exist on both sides and the local and macroscopic fields are related via the induced polariza- 
tion field [74] -see also [14] for a microscopic rigorous approach. In either case one obtains 
Etl = C,&J f , m the small cavity limit. In the case the emitter is a dipole equivalent to 
all the rest, the cavity is virtual and Cf = Cll = €M ° +2 is the Lorentz-Lorenz local field 
factor [55]. For the case that the cavity contains an emitter of much weaker polarizability 
than that of the host scatterers, the cavity is assumed to be empty and Cf — C emp = 
[56]. Same approach is followed in more complicated dielectric configurations [75]. 

It is obvious that the above procedure cannot be the correct one as it ignores the modes 
k > and hence the (^-dependent terms in \ J DOS^ niss ' lon . Nevertheless, let us assume 
for the moment that the above procedure is legitimum. The first erroneous assumption 
which appears implicitly in Eq.(210) is that ...the local field factor appears squared, as F (or 
LDOS) can be expressed in terms of and expectation value of the product of the two electric 
field operators (E Loc ) [76] - also in [22, 50, 77, 78]. The point being that the fact that the 
classical values of Ef-J and E e J* are related by a constant of proportion Cf, does not imply 
by any means that there is a constant of proportion C 2 between their quadratic fluctuations 
in vacuum, 

&dHf) = CM S {f) * (E^{r) ■ E?l f \r)) = C) (E e J f (r) ■ E e J"(r)). (211) 
The proof of Eq.(211) was given in Section VI B 3. We found there that the propagator of 



^ioc rea ds 



e// 

g$(k) = Um.{g l± }\ k ^ = ^Gf L {k), (212) 
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with ^ = C LL = while G e / f defined in Eq.(147) is the propagator of E e J f . Ne- 

glecting extinction, SjeMG 1 } <C ^{cmg}, and in application of the fluctuation-dissipation 
theorem, 

= £ LL (^(f) • E e J f \r)) oc £ LL n, (213) 

which is proportional to the MG coherent emission computed in Eq.(152) instead. An alter- 
native manner to understand where the confusion resides consists of interpreting the macro- 
scopic vector fields as operators acting on the coherent vacuum defined in Section VI B 1. In 
that framework, two factors Z 1 ^ 2 = C^l appear in the two-field vacuum expectation value 
in passing from (E e J f {f) ■ E e D ff \r)) to (E^J '(f) ■ E e ^p(r)). That is, 

(4 e //(r) • E\lp (f)) = Coh (Q\k ff {r) ■ E* eff \r)\tt) Coh = S - L (\Q\^ff (r) . | e ^t(^|Q)^- 
= C LL (k D ff (r)-k D m (r)). (214) 

The non-necessary implication expressed in Eq.(211) was first appreciated by the authors of 
[38] . The key argument given there being the lack of noise polarization. The authors of that 
work introduced an operator to account for such an effect and then required consistency. 
In our approach, we compute exactly the Green's function of the self-polarization field. 
Application of the fluctuation-dissipation theorem takes care of the polarization noise in an 
exact manner. 

In some old works it is claimed that it is only the bulk density of states, ~ n, the one 
responsible for the spontaneous decay which enters Fermi's Golden rule. That is based on the 
macroscopic quantization carried out for the first time in [79]. Such interpretation has been 
modified since then as further studies on the role and origin of the local field factors have been 
performed. However, still recent works [35, 49, 80] appeal to the erroneous interpretation 
that the density of states accessible to the photons emitted according to Fermi's Golden 
rule is just given by the bulk term. The latter corresponds to density of bulk normal modes 
instead which we have denoted by LT30S* ource/ess . The relation between the coherent (and 
only coherent!) emission spectrum and LDOS* ource/ess was made clear in Section VI B 2. 
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B. The actual nature of radiation and the need to go beyond the effective medium 
approximation 

It is a common error in the literature to identify the total radiative emission with trans- 
verse emission. As a matter of fact, the expression in Eq.(210) is intended as proportional 
to transverse spontaneous emission, eg. [19, 38, 81]. This error is related to the use of 
macroscopic fields and so to the neglect of short-wavelengths. 

In the first place there is a problem in imposing continuity conditions at the cavity sur- 
face. This is particularly problematic in the VC scenario. On the one hand, one assumes 
that an effective medium exists on both sides of the cavity, which implies some length scale 
over which the susceptibility is averaged in space. On the other hand, by imposing some 
'matching' conditions at the cavity surface one is assuming that the width of the surface 
is negligible in comparison to the length scale over which the spatial average has been per- 
formed on each side. However, that width cannot be less than the typical correlation length 
between the scatterers, £. That is precisely the radius of the inner cavity in which the spa- 
tial average is performed. Therefore, matching conditions are expected to fail in general but 
for those coherent modes which propagate in the 'effective' medium. Those are, the ones 
in Wf oh - and T^° h - on Eqs. (165,152). Because the rest of the transverse emission contains 
^-dependent terms, it is clear that Eq.(210) cannot be fully transverse. However, we have 
seen in Section VII A that Eq.(210) yields almost the right result for the ^-independent 
total decay rate up to order two in x -Eqs. (187, 188). This accidental coincidence is due the 
addition of the longitudinal dispersive terms of Eq.(183). We showed in the single-scattering 
model that, at leading order, the ^-independent longitudinal term equals the contribution 
of a local field factor. The resultant radiation was proved to be incoherent -see Section VI. 
Nevertheless, it is still possible to distinguish between ours and the usual formula without 
going to higher orders in x- If instead of using an integration sphere to collect the total 
radiation in far field the coherent field were measured, the difference could be verified either 
in experiments or in numerical simulations. 

One possible way to go around the problem with the boundary conditions would be to 
impose them not over the classical fields but over the Green's functions. In the VC model 
this is not a solution as one would face again the problem with the length scales. In the 
RC model things can be more promising provided R ^> £. However, one has still to assume 
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some form for the Green's functions on both sides as done in [82]. Unfortunately, in doing so 
all the intermediate correlations between the emitter and the medium are neglected. Only 
the first scattering after the emission of a virtual photon and the last scattering before its 
absorbtion get correlated to the emitter this way. 

It is also an error to attribute to the whole absorptive (^-dependent decay rate a longitu- 
dinal nature. This was already noticed in [82]. In the overlap approximation, some of the 
absorptive terms of Eq.(185) correspond to the transverse 7-factor of Eq.(171). 

C. Why the approach of [57] is not microscopic enough 

We comment on the approach of [57] specifically because it is referred by some authors as 
a microscopical proof of the usual formulae of the RC and VC models based on Eq.(210). In 
the following and for the sake of brevity we will not write all the equations in [57] but only 
those essential for our arguments. To avoid confusion, we will quote them within double 
brackets. 

The same as ours, the approach there bases on the computation of the Green function 
of a system which consists of an impurity placed within a host medium. They calculate the 
renormalized polarizability of the impurity following a similar procedure to that we used 
in Section III B. The host medium considered in [57] is a cubic lattice with lattice spacing 
£ m p -1 / 3 . The medium is therefore strongly correlated. However, because for the frequen- 
cies of interest k ^ <C 1, an effective medium can be defined with Xeff — Xmg for the reasons 
argued in Section VA. This is used by the authors to extrapolate their results to random 
media. The propagator of the local field is denoted by Q m in Eq.[(13)]. The authors first 
compute a transference matrix and a propagator denoted by Q, in the long- wavelength limit 
fc£ <C 1. The resultant propagator Q(R,R') of Eq.[(ll)] links any two distant lattice sites 
and contains only long-wavelength modes. Afterwards, the authors proceed to compute the 
t-matrices of impurities which are embedded in the host medium. To do so, they calculate 
Gm(r m ,r m ), where r m denotes the position vector of the impurity. In that calculation, the 
propagator of Eq. [(H)] is attached to the impurities in two different manners, depending 
on whether the impurity occupies a lattice site -in which case it is said substantial, with 
subscript s- or not -in which case it is said interstitial, with subscript %. That is done in 
Eqs. [(16,17)] through the computation of the 7^-matrices. 
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In the first place, the approach is not microscopical enough as the propagator in Eq.[(H)] 
does not contain the necessary resolution to probe distances shorter than £, and the distance 
from either an interstitial or substantial impurity to the first neighbors in the lattice is < £. 
Therefore, only transverse long-wavelength modes can be accurately accounted for. 

Nevertheless, for the case of an interstitial impurity their approach may be a good approx- 
imation under some conditions, as the impurity is connected through free space propagators 
to the T-matrix of the lattice and in the resultant function Qi(fi, fj), r*j is not a lattice vector. 
However, the correlation of the impurity to the lattice that way reduces to the correlation to 
the nearest lattice scatterers (the ones at the extremes of the T-matrix). An exact treatment 
would require the implementation of additional correlations as we did in Section VB. 

The situation is worse for the case of a substantial impurity. This can be seen by com- 
parison with the exact result we obtained in the virtual cavity scenario. That is, for the 
case that the polarizability of the substantial "impurity" is equivalent to that of the rest 
of particles of the lattice, the resultant setup would correspond to the strict virtual cavity 
scenario. Therefore, in the long-wave length limit, one would expect to obtain a result pro- 
portional to Eq.(152) for the renormalized long- wavelength decay rate -which is identified 
in [57] with the radiative decay rate after Eq. [(15)]. Instead, they obtain an additional local 
field factor which supports the erroneous introduction of one local field factor per electric 
field operator. The reason for the mistake is that the authors make an erroneous usage of the 

— * — * 

T-matrix and the propagator Q(R,R') of Eq.[(ll)]. The T-matrix in Eq.[(6)] just includes 
terms which link any two dipoles in the lattice. However, if one is to compute G(R, R) out 
of it, self-correlations must be explicitly considered. That is, negative correlations are cor- 
rectly implemented in Eq.[(6)] by excluding the vector R = in all the summations, which 
leads to the equivalence with the overlap approximation in a random medium. However, 
when distant scatterers are also correlated -eg. through self-correlations- such an addi- 
tional correlation must be incorporated as conditions over the lattice vectors of the sums. 
The translation invariance assumed in the sums gets 'virtually' broken when introducing 
self-correlation. That is the information which is missing in the approach of [57]. As a 
consequence, the evaluation of Q(R,R!) at R' — R does not incorporate self-correlation in 
the right manner and its use in Eqs. [(14,16)] is incorrect. 

In the following, we analyze the formulae of [57], compare them with ours and spot the 
errors. We have illustrated the failure in Fig. 19. Eq.[(6)] with the approximation of Eq.[(7)] 
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FIG. 19: (a) Diagrammatic representation of the T-matrix as computed in [57], Eq.[(6)], in the 
long- wavelength limit of Eq.[(7)], 1. That limit is effectively equivalent to taking the overlap 

approximation, (b) Diagrammatic representation of two 'apparently' different diagrams in Q(r,r). 
They come to be equivalent when self-correlations are implemented. In the overlap approximation 
of (a), because the 1PI piece underlined on the r.h.s. is not included in the effective self-energy 
function under the overlap approximation, such an equivalence disappears in the approach of [57]. 
(c) Feynman rules. 

is equivalent to 

T(k^O) = -k 2 a\l + k 2 pa[G {0) {k) + C{k£^0)]\~\ (215) 

where C(k£ — > 0) is the cavity tensor with components given in Eqs. (90,91). In the long- 
wavelength limit, C(k^ — > 0) = — 3^7 ^- a — a o[l + ^^ a o] i s the in-free-space radiatively 
corrected single particle polarizability -the authors of [57] denote it by a instead. Only 
the two-point correlation function, hc(r) = — 0(r — £) with £ p -1 ^ 3 , enters the cavity 
factors. That implies that, effectively, the same negative correlation function which enters 
the computation of Xmg m a random medium does so for the computation of the T-matrix in 
a cubic lattice in the limit k£ <^ 1. Expanding next the denominator of T as a perturbative 
series in k>lpa[G^ (k) + C(k£ — > 0)], it is easy to verify that the overlap approximation 
is implicit in this expansion such that the susceptibility there is that of MG in the limit 
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ki -> 0. Thus, f eff (k) = reads, 

f eff (k) = p-'ixMG + XM G G eff (k) X M G }, (216) 

where G e ff(k) is given in Eq.(147) and is denoted by Qdiei m [57]. Hence, for any two distant 
lattice points R ^ R', Eq.[(ll)] can be otained in the long-wavelength limit by removing the 
bare vertices a , pa attached to the extremes of Eq.(216), 

^ = ( a / a o) 2 C 2 LL g diel (R, R'). (217) 

If one is to compute Q e ^ (R, R) using the general expression Eq. [(4)], the T-matrix there is 
not that in Eq.[(6)], which corresponds to the one which enters the Dyson equation Eq.(56) 
for a cubic lattice. The point being that that in Eq.[(6)] was computed without considering 
the possibility that the first and the last dipoles in the diagrams were the same. As a matter 
of fact, application of the overlap approximation from the very beginning prevents from 
correlating two distant scattering events -see Fig. 19- and hence the diagrammatic 'trick' we 
used in Section IV B to deduce the exact formulae turns inapplicable. Hence, if one aims to 
write the equation for Q(R, R) in the form of Eq.[(4)], one should use a modified T-matrix 
whose form derives from our Eq.(63), 

Tui(k) = -[G^ril + ^T UI ]. (218) 

Because a lattice is a strong correlated system, the exact formulation of the problem in 
terms of the t-matrix is more convenient as the correlation length can be extended eventually 
further than £. The formula to compute the propagator Q(R, R) is that in Eq.(73), 

g(R, R) = zlj2 G (0) (n) • m, (219) 

Note also that, for the sake of consistency, the fully renormalized a polarizabilities enter 
this formula, and neither the in-free-space radiatively corrected one, a, nor the static one, 
a - The components of of Eq.[(6)] in [57] relate to t(fj) in our nomenclature through 

I + fcjja^G^rOe**] -1 , {f<} = £(ai + bj + ck), a, 6, c G Z, n ^ 0. 

(220) 



f s = -k 2 a 



{*} 
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Writing the t-matrix in the alternative form, T(f) = ^{rj t(fi)8^(f — fi), Eq.(219) reads 

S( R, A)^l dV ■ m = g / |* G<°> (*) • f, (22!) 

For the long- wavelength modes, the same result as in an MG dielectric holds at leading order 
in £; £, 

G(R,R)\*«i = ^1 ^^v^- (222) 

07T O 

For higher frequency modes, it is unavoidable to compute the first terms of Tg. 

X. THE ZERO-TEMPERATURE ELECTROMAGNETIC PRESSURE 

The role of the zero-point-energy (ZPE) in a variety of physical phenomena has received 
renewal attention in the last years. In particular, a strong motivation is to estimate its 
contribution to the cosmological constant, A, in connection with the possibility of being 
measured in the laboratory. The need for the existence of the cosmological constant is itself 
a controversial matter which is out of the scope of this paper -see eg. [83-85]. Standard 
cosmology predicts a small value for the current A, which needs of an apparently unnatural 
fine-tuning. Leaving aside some exotic contributions to dark energy, at least two main 
contributions to A must be considered. The first one comes from the Standard Model 
of particles in the form of condensates of quarks and gluons and the Higgs itself. The 
second one comes from the ZPE of quantum fields. Eventually, these energies give rise to 
negative pressure which sources the repulsive gravity force that leads the expansion of the 
universe. The main difficulty in dealing with ZPE computations is that the integrals involve 
ultraviolate divergences. One could think of the need of being regulated by some physical cut- 
off. However, this is not yet fully justified [86, 87]. Nevertheless, it is customary -as quoted 
in [2]- to present genuine quantum electromagnetic effects like the Casimir effect, the van der 
Waals forces, the Lamb-shift, etc. as evidences of the existence of the ZPE fluctuations of the 
EM field. Further on, some papers have suggested the possibility of measuring the ZPE of the 
electromagnetic field in Josephson junctions [88]. As pointed out by Jetzer and Straumann 
[89, 90], the only physical observable which can be measured in absence of coupling to gravity 
is the variation on the vacuum energy as a response to external couplings. Any physical 
observable measured this way remains invariable to any prescription on the normalization 
or regularization of the vacuum energy. Hence, the scheme we followed to regularize the 
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shifted-electrostatic polarizability, oP stat , was merely phenomeno logical. Infinite quantities 
were swept under the carpet and remained unaffected by the rest of the renormalization 
procedure. That regularization scheme only served to us to accommodate in the same 
footing further finite contributions. It is analogous to the regularization of the divergences 
in the electron self-energy which appear in QED in computing the experimental mass of the 
electron [29, 30]. Figs.5(a, b, C2, ofe) depicts how, after formal integration of those EM modes 
which give rise to the bound atomic state A, one can interpret that the 'remaining' vacuum 
fluctuations which amount to the electromagnetic ZPE interact with such a state giving rise 
to radiative corrections. However, this is just a result of the order in which fluctuations are 
integrated out. It is plain that there is no transfer of energy from the ZPE into the dipoles 
self-energy. 

Jaffe pointed out in [2] that the Casimir force can be computed out of the derivation of the 
density of electromagnetic states with respect to some parameter on which the background of 
matter fields depend. In the case of the force between two metallic plates, the background 
fields are the electric currents confined on the plates separated by some distance. The 
diagrams involved contain external legs representing the fields on the plates -those are 
the current loops depicted on the left plate of Fig.20(c). This way, the diagrams with 
photons attached to the plates cannot be interpreted as vacuum-diagrams. An alternative 
interpretation is given by Milton et al. in [91]. The authors identify the interaction energy 
between the plates with the self-energy which renormalizes their (divergent) bare masses. In 
Fig. 20(c) the photon propagators carry the radiative corrections which give rise to the plates 
self-energy. In the interpretation of [2], it is the variation of the density of those photon 
states which yield the Casimir force. In addition, in the context of Lifschitz formalism [92] 
Milton identifies [1, 93] a bulk energy which depends on the dielectric medium between the 
plates and is proportional to the volume. 

We proceed next to describe where energies and forces reside in our formalism. It was 
emphasized throughout Section II that actual dipoles play a double role as polarizing the 
EM vacuum and renormalizing their own polarizability. In energetic terms, this gives rise 
to a change on the energy carried out by the fluctuations of sourceless EM modes and to a 
shift on the resonance frequency of the dipoles. The former is the analog of the bulk energy 
of [93] while the latter is the analog of the self-energy of the metallic plates of [91]. We will 
denote the energy density on sourceless EM modes by T ra d and that stored in the dipoles 



82 



by Tmat- The corresponding bare energy densities are the divergent ZPE of Eq.(lO) and 
To = phuA for dipoles in the atomic state A. While the latter, associated to actual matter, 
can be detected in EM phenomena, the former eventually needs of coupling to gravity to be 
detected [94]. 

Once the dipoles are considered bound states, statistical translation invariance holds 
at length scales greater than a. It is in this sense that we can talk of a proper vacuum, 
in contrast to the scenario of the parallel metallic plates. As it was illustrated in 
Section II B, the diagrams contributing to the energy of the polarized EM vacuum are 
topologically identical to those contributing to the dipole self-energy. The only difference 
being a change on the reference frame. Fig.20(a) shows a typical diagram of both T ra d 
and Tmat- It contributes to T ra d when all the points on the diagram are equivalent. It 
contributes to T ma t when the origin and end of the diagram is chosen to be a point scatterer 
as in Fig. 20 (6). The main difference between T ra d and T mat is that while the former can be 
computed at zero-order in perturbation theory together with JF , the remaining of T ma t is a 
second-order contribution. That is, 

N „ N 

T rad + T = (n^A(r) \V~' / d\ e j\E(f)\ 2 + V~ 1 H (f)\j2^A(r)^) (223) 

3=1 J J=l 

= |^(fi|||(f)| 2 |fi)^-+^ = -^A_ J c 2 Tr{3{G ± (f,f)}}dc + ^ (224) 

= c' 1 J fvuj LDOS* ourceiess dcu + T , (225) 

where H is given in Eq.(ll) and V is the volume occupied by the dielectric and N is 
the number of dipoles. The atomic wave functions of the dipoles in state A have been 
explicitly separated so that their contributions to To are additive. The radiative term on 
Eq. (224) assumes knowledge of x± f° r an y frequency. However, we have already argued that 
our approach has a natural frequency cut-off at c/a. Also, in general, the dipoles contain a 
number of atomic resonances so that To should contain a sum over those resonances instead. 
For the sake of simplicity and in order to give a closed formula, we will restrict ourselves 
to the computation of T ra d up to frequencies where an effective medium exists, u < c/£, 
yielding, 

FradUc/t + J dcu LU 3 [n(Lu) - 1], (226) 

where n{uj) = 3?{^/e^}. LDOS in Eq.(225) is referred to in some papers as radiative density 
of states [35, 49]. However, as stressed in the previous section, it must not be confused with 
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that entering Fermi's Golden rule. The second term in Eq.(226) is only non-zero within the 
volume occupied by the dielectric. Hereafter we will denote it by Tf^d = Frad — £°- 

For the computation of the remaining T mat we have to go to second-order in perturbation 
theory to calculate the real part of the self-energy of each dipole. Formally, it reads 

» N N 

T ma t-T* = V~ l / d 3 r J2 ^{(n,J2^A(r) \H int (f,t)\ ^^(r),7^) (227) 

N N N 

x (ZX(*D,7«* \Hint(f,t)\ ^Vi(r1,^)n[^-^-M _1 } (228) 

j=l j=l 1=1 

I,U! 

x [hu A - tvjj! - (229) 



where Hi nt is given by Eq.(13). We will simplify matters by considering just two-level dipoles 
so that Ij in the previous equations can take only the values A or B, with Uq = uja — <^b 
the resonance frequency and /x the transition dipole amplitude. Next, we make use of the 
renormalization scheme developed in Section III B and identify 

Tmat = phc Ak res = ^M 2 P9fc{2 7 * + i\}\~ k=kres . (230) 

Ak res and k res deduce from Eq.(47). Our renormalization procedure takes care of possible 
divergences in the computation of Eq.(229). Hereafter, we will denote the dielectric self- 
energy by T s ^ a { E = Tmat — At leading order in pa , the van-der-Waals energy J^ vdW 
between two random dipoles in an MG gas is given by Eq.(230) replacing 9ft{27^ + ^}\ k=kres 
with the leading order term of ^{^{k)} given by Eqs.(116,117). That energy corresponds 
to the electrostatic interaction between two fluctuating random dipoles, 

r dW = !^|2pft{ 7 |M)} 



3e 



o 

= P^H^o Pie = "P ^"oP/e (231) 

The mutual induction is given by the 2 nd and 3 rd diagrams of Fig. 14(6) with the propagators 
carrying only longitudinal photons. 

The dielectric and the corresponding polarized vacuum are characterized in our case by 
the density of dipoles and a typical correlation length £ between them. The former plays an 
analogous role to that of the separation between the metallic plates in the Casimir setup. 
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In a van-der-Waals gas, £ is the van-der-Waals radius. Therefore, internal forces between 
dipoles originate as a response to variations in both/either p and/or £. In high correlated 
systems they both are related. For a gas, £ is fixed and the pressure reads, 

d 

Ptot = ~gy{V C^Vod + Fmat )} 

a -pdiel a <pselfE 

= [P^jf- - + iP^p- - F™t E \ = Vrad + V mat . (232) 

Because the medium is statistically isotropic and homogeneous, the pressure is a scalar which 
tends to expand or contract isotropically the distance between particles. Note that only V ma t 
is considered when introducing van-der-Waals-Casimir forces in the free energy of a complex 
medium. No reference is made to the electromagnetic pressure, V ra d- Prad is interpreted as 
the pressure exerted over the dipoles by the vacuum fluctuations of the normal EM modes. 
The fact that T, Ak res and P ma t are all proportional to 7-factors and their derivatives implies 
that none of them can be attributed to transference of energy from EM vacuum fluctuations. 
They are rather associated to changes in the self-energy of material degrees of freedom. 

It is now straightforward to compute the zero temperature pressure of a van-der-Waals 
gas. That is usually written as [20] P^^ = —a'p 2 , where a' parametrizes the dipole-dipole 
interaction. Applying Eq.(232) to Eq.(231), 

6tt £ 3 

we identify a! = ^ao/£ 3 - Should we introduce the higher order terms in p of Eq.(179) we 
would be performing the virial expansion of restricted to van-der-Waals short-ranged 

forces going like ~ 1/r 7 . 

For the sake of completeness we compute in the same approximation V ra d- We will restrict 
ourselves to the simplified case in which the dielectric possesses a unique resonance with not 
too much loss of generality. Making use of Eq.(226) together with Eq.(195), we obtain 

Fral - ^2^oP (r + Ar)(A; + AA; res + Afc L ), (234) 

where only finite p-dependent contributions have been considered in the range of validity of 
Eq.(195). Next, using the values of the 7W-factors computed in Section VIA, 

ArW ~ ? -pa o r , AfcN, ~ ~^a 2 p/e, Ak L = -h o a p, (235) 

O D7T 3 

we obtain, J** ~ ^k 7 a 2 p [1 + pa (^ - (236) 
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V vdW = -^t|p 2 , (233) 



Applying Eq.(232), we get at leading order in pa , 

For realistic values of £ and cko, ^rad is positive except for extremely long- wavelength res- 
onances of the order of 10 6 nm. This means that the contribution of the individual self- 
energies of the dipoles in V^} d is negligible in comparison to that coming from the additional 
decay rate, ArW, and the Lorentz-shift, Ak L . The latter being just proper of the effective 
medium. Because A/c L enters with opposite sign w.r.t. that of Ar^, its net effect is to 
subtract from Ar' 1 ' the longitudinal contribution of a local field factor at leading order in 
pa , ArW + Ak L = ArJJ 1 . The ratio between Eq.(237) and Eq.(233) is, 

V W _ 5 

l^ri 10 2 . (238) 



<pvdw i28n 2 

for realistic values of a . Therefore, we find that V ra d is negligible for a van-der-Waals gas 
in comparison to V mat . 

We finalize this Section speculating about possible observational effects of the above forces 
at cosmological scales. It is known that the magnetic dipole associated to rotating galaxies 
is strong. Therefore, the interactions between distant galaxies or clusters of these which are 
not gravitationally bounded might be predominantly electromagnetic. As long as galaxies 
can be treated as point dipoles, we can model the universe as a cold van-der-Waals gas. Its 
pressure, disregarding the coupling to Hubble's flow, would be determined by the van-der- 
Waals-Casimir forces and the vacuum radiative pressure. However, we must emphasize that 
this does not imply by any means that the associated energy T s ^ t E + must be of any 
gravitational relevance. Said this, it results tentative to associate the negative pressure V^j&Q 
-which however only in flat space and ignoring retardation effects has been proved here to 
be negative- with a cosmological constant. Nevertheless, V to t might generate a red/blue- 
shift of similar characteristics to the cosmological one at scales where the distribution of 
galaxies is really isotropic and homogeneous. Otherwise, it would yield small variations on 
top of the cosmological red-shift, due to the generation of peculiar velocities. Evidently, 
curvature, causality, magnetic effects and retardation effects, so far ignored in the present 
work, should be considered. Also, we would like to stress that, although we have showed 
that ordinary QED effects are explained without reference at all to the ZPE, this does 
not imply that current investigations on more exotic Casimir-like effects are meaningless. 
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Note that our derivation of the total pressure in Eq.(232) considers only variations on the 
configuration of actual dipoles, which does not affect to either JF or £°. However, other 
kinds of variations may be thought of which so do. Hence, in the so-called Casimir effect in 
compact dimensions and the dynamical Casimir effect in a manifold with moving boundaries, 
the resultant pressure is associated to variations on the geometry of the space-time manifold. 
That gives rise to the aperture/clousure of channels for the normal modes of quantum fields 
and hence to variations on £°. The fact that the 'ordinary' Casimir effect between parallel 
plates can be explained effectively in similar terms is just due to the fact that the interaction 
between perfectly metallic plates can be well approximated by boundary conditions over 
the mediating EM field which 'mimic' the effect of a compact space [2]. However, those 
conditions are imposed over Q and hence give rise to changes on the self-energy of the plates 
[91] and not over £°. In this sense, V ma t for a gas of galaxies is the actual cosmological 
analog to the ordinary Casimir's force. 



XI. CONCLUSIONS 



We have studied the phenomenon of electric dipole emission in statistically homogeneous 
complex media. Based on the computation of the spectrum of sourceless normal modes and 
the spectrum of the EM fluctuations which yield the self-energy of the dielectric matter, 
we postulate the distinction between a sourceless vacuum, |fi) sX , and a self-polarization 
vacuum, \Vt) s ' p \ The former is a proper vacuum as it is translation- invariant. The latter 
is defined upon each dipole. In physical grounds, the distinction bases on the observation 
that the constituents of a dielectric play a double role. In the first place, they polarize the 
zero-point vacuum, |0). Complementarily, they renormalize their own polarizability, cto, for 
the multiple-dipole interactions generate a self-energy. The EM energy in \ fl) s - L is computed 
non-perturbatively. The dielectric self-energy enters at second-order of perturbation theory 
following Fermi's Golden rule. 

In application of the fluctuation-dissipation theorem, the spectrum of EM fluctuations of 
\Q,) S - L is determined by he Dyson propagator, G. The spectrum of fluctuations of the self- 
polarization field in \Q) S ' P ' is determined by Q. The latter satisfies a Lippmann-Schwinger 
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FIG. 20: (a) Typical 1PI diagram contributing to both Tf^ and T s ^ t E . It presents a characteristic 
fractal structure, (b) Diagrammatic representation of LT>OS < ^ nisslon . Translation invariance around 
the loop is broken at the position of the emitter, (c) Sketch of the Casimir self-energy which 
renormalizes the mass of two perfect conducting plates. On the l.h.s., current-loops are depicted. 
On the r.h.s. the representation mimics the formalism we used for the radiative corrections which 
renormalize the bare polarizability of small dipoles. 

equation with a stochastic kernel. In a homogeneous dielectric the kernel is given by 



3ISW = 



pa 



X±,\\(k) G®,(k) 



Xi#) 
pa 



As a result, 



X±(k) 



pa 
1 



pa k 2 [l + x±(k)}-k 2 ' 



1 



gr(k) = _ xll (fc)Gn(*) = — 



X|| (k) 



pa " ' pa fco[l + X||(*0]' 

The above equation is exact in the strictly virtual cavity scenario (VC). For the self- 
polarization of an impurity, the computation of Q is model-dependent. The model is gener- 



ically referred to as real-cavity model (RC). It is possible to compute the kernel of the 
equation for Q in closed form for the case the impurity is placed in a large cavity with 

K ±,\\( k ) X±,\\(k) GY :ll (k) 1 

where n^ik) is derived formally in Section VB. For a medium in which dipoles remain 
fixed Q can be expressed in function of the transference matrix as, 



Q(r , f ) = — ^ ^ G(°)(r , r,) • i(f h f ). 



The local density of EM states in is denoted by L£)OS^ m * sslon for it corresponds 

to the density of channels accessible to the emitted photons. The power emission, W u , the 
decay rate, T and LDOS^"™ SSSMm are decomposed according to their transverse/longitudinal 
and coherent /incoherent nature as, 

2W^ h - = W Q [ ^2^{^}-s{G ± (k)}, 
J (27r) d pa 

2W Incoh. = w f ^_ 2 Q { X±i^_ mG±{k)h 

J (zny pa 

W{ nc ° h - = W [ ^3{^}»{G,|(*)}, 
11 J {Ztt) 6 pa 

where W a must be substituted by the appropriate constants in the case of T and 
LDOS^ misss * on . We found that the vacuum in which coherent emission propagates, \Q) Coh , is 
additionally polarized w.r.t. that of bulk normal modes, |f2) sJ \ The renormalization func- 
tion being, Z± t \\ = In a Maxwell-Garnett dielectric, it equals a local field factor. 

We have computed the decay rate of a weakly-polarizable interstitial atom in an MG 
dielectric. The results are meant to fit the experimental data on the life-time of Ce +3 ions. 
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We obtained, 

1 MG — 

r \UDis P er. = 2 f 3?{i XMG + l - x 2 MG - 0.051 Xmg + 0-055 XmG ~ 0-015 Xmg + •••} 
2f 1 1 

+ ^{ 2^mg + ^Xmg ~ 0-076 x^g + 0-083 X mg ~ 0-022 X mg + •••} 

2f 11 

+ ^3{Xmg + gXl/G - gXL? + 0.073 x 4 MG - 0.028 XmG + •••}, 

where r o contains the renormalization of //. In the single-scattering approximation, a clas- 
sification of the decay rate in terms of radiative and non-radiative energy transfer has been 
given -see Fig. 16. 

The computation of the dielectric constant of an MG dielectric is a problem of self- 
consistency due the double polarization role played by point dipoles. 

We have computed the decay rate of a low-polarizable substantial impurity within a 
large cavity. The results are meant to fit the experimental data on the life-time of Eu +3 
ions. Formulae are exact up to order C(x 2 ), 

r c R G h - = r [K{^{(^) 2 }-^{ e -i}\ 
rK— = w{ e ~ {e ~ 2)2 }> 
r$T 6 ' = - r oHv^{(^) 2 } 

/l 1 \^ r 2(e-l) 2 , ,,i 

In comparison to previous works, we have found that the usual formulae for the decay 
rate in the virtual and the real cavity models are erroneous in physical grounds. The main 
reason being that the inclusion of one local field factor per field operator entering Fermi's 
Golden rule is an erroneous assumption. Also, only coherent emission can be accurately 
computed using macroscopic fields and incoherent radiative emission is longitudinal. 

The total EM energy density is the sum of the energy of the EM fluctuations in |fi) sX , 
Frad — I ^ 2 Tr j (f, f) } jdc<j, plus the total self-energy stored in the dielectric, 

Fmat — Fo + ^|/ i | 2 ^ 2 3^{27i + 7|| }|fc =A , res . It is the variation of the latter w.r.t. external 
parameters, £ and p, that gives rise to generalized van-der-Waals-Casimir forces in a 
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complex media. Neither zero-point EM vacuum energy nor JF contribute. The variation 
of T r ad + T m at w.r.t p yields an isotropic pressure. The zero-temperature pressure of a 
van-der-Waals gas is found to be J> vdW = — ^ck 2 , p 2 /£ 3 - Higher order terms of the virial 
expansion have been found together with additional contributions from long-ranged forces. 
The radiative pressure at zero temperature is estimated to be at least twelve orders of 
magnitude weaker at leading order in p. 

We have argued on the possibility of testing our formulae for e e ff, T and resonance 
shifts in experimental setups. Cold atoms and colloidal liquids are suggested to be suitable 
candidates. In order to distinguish ours from the usual formulae of the virtual cavity and 
the real cavity models, it has been argued that either high refractive index or optically thick 
media are needed. Alternatively, measuring the coherent radiation instead of integrating 
the total one may do. 

The detection of the zero-temperature radiative pressure, which is estimated inviable 
in gasses, might be possible in highly correlated cold fluids which present structural 
resonances. However, we need a more precise dielectric function than that for the two-level 
mono- atomic dielectric used in the derivation of Eq.(237). The dielectric constant should 
be valid for a much wider range of frequencies. To this respect, Kramers- Kronig relations 
[95] might help in performing the integration of Tf % ad l . More investigation is needed on this 
point. Finally, predictions at cosmological scales need of the introduction in the present 
formalism of curvature, retardation and magnetic effects. Work on this matter is in progress. 
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NMP-2007-SMALL-1 and the Scholarships Program 'Ciencias de la Naturaleza' of Ramon 
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